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ABSTRACT 

We present the color-magnitude and color-stellar mass diagrams for galaxies with Zphot ^ 2, based 
on a i^(^B) < 22 catalog of the \ x \ U° Extended Chandra Deep Field South (ECDFS) from the 
MUltiwavelength Survey by Yale-Chile (MUSYC). Our main sample of 7840 galaxies contains 1297 
M« > 10^^ M© galaxies in the range 0.2 < Zphot < 1-8. We show empirically that this catalog is 
approximately complete for M, > 10^^ M© galaxies for Zphot < 1-8. For this mass-limited sample, 
we show that the locus of the red sequence color-stellar mass relation evolves as A(u — r) cx (— 0.44± 
0.02) Zphot for -Zphot ^ 1-2. For Zphot ^ 1-3, however, we are no longer able to reliably distinguish red 
and blue subpopulations from the observed color distribution; we show that this would require much 
deeper near infrared (NIR) data. At 1.5 < Zphot < 1-8, the comoving number density of Af* > 10^^ 
M© galaxies is ^ 50% of the local value, with a red fraction of « 33 %. Making a parametric fit to the 
observed evolution, we find ntot{z) oc (1 -I- Zphot)""'^^^"'^^^^"'^"^- We find stronger evolution in the 
red fraction: frod{z) oc (1 -I- Zphot)^^'^''^^'^^^''^'''^"^^- Through a series of sensitivity analyses, we show 
that the most important sources of systematic error are: 1. systematic differences in the analysis of 
the z and z ^ samples; 2. systematic effects associated with details of the photometric redshift 
calculation; and 3. uncertainties in the photometric calibration. With this in mind, we show that our 
results based on photometric redshifts arc consistent with a completely independent analysis which 
does not require redshift information for individual galaxies. Our results suggest that, at most, 1/5 
of local red sequence galaxies with > 10^^ Mq were already in place at z ~ 2. 
Subject headings: Galaxies: Formation — Galaxies: Evolution — Galaxies 
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1. INTRODUCTION 

Observing the evolution of the massive galaxy popu- 
lation provides basic constraints on cosmological mod- 
els of structure formation, and so helps to identify the 
physical processes that govern the formation and evo- 
lution of massive galaxies. In this context, the color- 
magnitude diagram (CMD) — astronomy's most basic di- 
agnostic plot — has been particularly important and use- 
ful over the past five years. Physically, a galaxy's rest- 
frame color is determined by the (luminosity weighted) 
mean stellar age, modulo the mean stellar metallicity 
and extinction from dust in the ISM. The restframe op- 
tical brightness acts as a proxy for the total stellar mass, 
although the connection between the two has a simi- 
larly complicated dependence on star formation history, 
metallicity, and dust. The CMD thus offers two com- 
plementary means of characterizing the star formation 
history of individual galaxies, in terms of the amount 
and character of their starlight. 

In the local universe, galaxies can be separated into two 
distinc t but overlapping p opulations in color-magnitude 
space ()Baldrv et al.ll2004 ): a relatively narrow and well- 
defined 'red sequence', as distinct from the more diffuse 
'blue cloud', with each following its own color-magnitude 
relation (CMR). Red sequence galaxies dominate the 



bright galaxy population, and tend to have the more con- 
centrated light distr ibutions typical of morphologically 
early type galaxies ("Stratev a et al.l 120011: iBlanton et all 
[20M lDriveriraI] r200f): Va n der Well 120081 ). Thev tvp- 
ically have stellar masses greater than 10^°'^ M© and 
are dominated by old stars, whereas blue cloud galax- 
ies are typical ly less massive and con t inue to be actively 
star forming (|K"auffmann et al.l [20031 : iBrinchmann et aTl 
120041: IWvder et al.ll2007f) . Further, red sequence galax- 
ies lie pref e rentia l ly in higher density environments 
(iHogg et all 120031 fBlanton et all l2005at iBaldrv et all 
l2006l ) . The emergent picture is of a population of mas- 
sive, quiescent, concentrated, and strongly clustered red 
sequence galaxies, as distinct from the typically less mas- 
sive, disk-dominated, and star forming blue cloud popu- 
lation (Ellis et al. 2005). This paper focuses on the red- 
shift evolution of the red sequence galaxy population. 

Using high qualit y photome t ric red shifts from the 
COMBO-17 survey. Bell et al.l (|2004bf ) showed that a 
red galaxy sequen ce is already in place at Zphot 
1 (see also, e.g., E ^eLaU 120021 : IWeiner et all l2005l : 
IWillmer et al.l 120061 ) . Further, as in the local universe, 
the Zphot ~ 0.7 red sequence is do minated by passive , 
morphologically early type galaxies l|Bell et al.l l2004aD . 
The combined mass of red sequence galaxies at 2 ~ 1 
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is at l east half of the present day value (|Bell et al.l 
I2004U iFaber et all l2007t iBrown et al.]l2008l ). By con- 
trast, the stellar mass density of actively star form- 
ing blue cl oud galaxies remains more or less co nstant 
for z < 1 (|Arnouts et all l2007l : iBell et all [2001 . even 
as the combined star formation rate drops by an or- 
der of magnitude over the sanie int erval (|Lillv et al.1 
ll996HMadau et al.lll996l : lHopkinsll2004D . These results— 
a steadily growing number of passively evolving red 
galaxies, and a relatively constant number of actively 
star forming blue galaxies — have led to the idea of 
a quenching mechanism for star formation, operating 
to incite a transformation that moves active galax- 
ies from the blue c loud onto the passiv e red sequence 
jMenci ct al 2005; Croton et al. 2006; Cattane o eFal 
[2006; Dckcl & Birnboim 2006; Dc Lucia ct al. 2007fK 

Our specific goal in this paper is to quantify the evolu- 
tion of massive galaxies in general, and of red sequence 
galaxies in particular, in the color-magnitude and color- 
stellar mass planes for Zphot ^ 2. The i ^ z <2 interval 
is particularly interesting: whereas the z 1 galaxy pop- 
ulation appears qualitatively similar to the local universe, 
at least in terms of the existence and properties of red 
sequence galaxies, the situation at 2; > 2 may be quite 
different. While massive, passive galaxies have be en con- 
firme d at z > 1.5 (iDaddi ct al, 2005; McGrat h et all 
I2007f) and even z>2 (|Kriek et al.ll2006f) . these galaxies 
do not appear to dominate the massive galaxy popula- 
tion as they do at z < 1. Indeed, it appears that the 
median massive galaxy at z ^ 2 has the infrar ed lumi- 
nosity of a LIRG or ULIRG (iReddv et al.l[2006h . More- 
over, whereas the number density of massive galaxies at 
z ~ 1 is > 50 % of the local val ue (lJu neau et al.ll200"5l : 
iBorch et al.ll2006l : iScarlata et al.r 2007). a t z > 2 it is in- 
ferred to be < 15 % fFontana et al. 2006; 'Arnouts et al.' 
I2007t iPozzetti et"a l. 2007; Percz-Gonzalcz et al 2008). 
This marks the redshift interval 1 < z < 2 as po- 
tentially being an era of transition in the universe, in 
which massive galaxies first begin both to appear in 
large numbers, and to take on the appearance of their 
local antecedents. This coincides with end of the pe- 
riod of peak star formation in the universe; while the 
cosmic star formation rate rises sharply for z < 1, it 
appears to pl a teau or even peak for z > 2 (see, e.g. , 
Hopkins! 120 04': 'Nagmine et al.' '2006^; 'Panter etafl l2007t 
Tresse et a l. 2007; Pcrcz-Gonzalcz et al. 200|). What- 



ever the mechanism that quenches star formation in mas- 
sive galaxies may be, it is in operation at 1 < z < 2. 

The technical key to gaining access to the 1 < 
z < 2 universe is deep near infrared (NIR) data 
(jConnoUv et al.|[T997D . since at these redshifts the rest- 
frame optical features on which both spectroscopic and 
photometric redshift and stellar mass determinations rely 
are shifted beyond the observer's optical window. More- 
over, the inclusion of NIR data makes it possible to con- 
struct steU^i^jMSsJirnite with high complete- 
ness ()Van Dokkum et al.l 120061 ). Among the next gener- 
ation of NIR-selected cosmological field galaxy surveys, 
the MUltiwavelength Survey by Yale-Chile (MUSYC; 
Gawiser et al. 2006) is among the first to become pub- 
lic. MUSYC has targeted four widely dispersed South- 
ern fields, covering a total of one square degree in the 
U BVRIz' bands. Coupled with this optical imaging pro- 



gramme, MUSYC also has two NIR imaging campaigns: 
a wid e {K'^^^^ < 22) component over thr ee of the four 
fields (|Blanc et al.ll2008l : iTavlor et al.ll2008l hereafter Pa- 
per I), and a de eper (K^-^^^ < 23.5 ) component for four 
10' X 10' fields (|Quadri et alll200l . 
This present paper focuses on the Extended Chandra 



Deep Field South (ECDFS), one of the four 



MUSYC fields. Centered on the historical Chandra Deep 
Field-South (a = 03 '^32'"28s, 6 = -27°48'30"; J2000 — 
iGiacconi et al.|[200ll ) , this is one the best studied fields on 
the sky, with observations spanning the full electromag- 
netic spectrum from the X-ray to the radio. N otably, this 
field is also a part of the COMBO- 17 survey (|Wolf et all 
l2003f ). and has received Hubble Space Telescope ACS coy- 
erage as part of the GEMS project (|Rix et al.ll2004D . as 
well as extremely deep Spitzer Space Telescope imaging 
from the SIMPLE (iDamen et al.ll2008D proje ct. Further, 
the GOODS project (jOickinson et al1l2002 f) covers the 
160 □' at the centre of this field, including support- 
ing NIR data from the ISAAC instrum ent on the VLT 
(|Grazian et al.ll2006t IWuvts et al.l 120081) . Complement- 
ing these and other imaging surveys, a wealth of spec- 
troscopic redshifts are available from large campaigns in- 
cluding the K20 survey (Cimatti ct al. 2002), the VVDS 
project (|Le Fevre et al.ri2004D. the two GOODS spec- 
trosc opic campaigns (jVanzella et al.ll2008l:[Popesso et al.l 
l2008l ). and the IMAGES survey (|Ravikumar et al.ll2006D . 
among others. 

The plan of this paper is as follows. We begin in 3Hby 
giving a brief overview of the data used in this paper — 
the MUSYC ECDFS datas et itself, as well as the z « 
comparison sample from Blanton et al.l (|2005bf) . Next, in 
^ we describe our basic methods for deriving redshifts 
and restframe properties for z ^ galaxies; our analysis 
of the z « comparison sample is described separately in 
Appendix [XI In fJU we construct a stellar-mass limited 
sample of massive galaxies from our if -selected catalog. 

Our basic results — the color-magnitude and color- 
stellar mass diagrams for Zphot < 2 — are presented in 
iJ51 We then analyse three separate aspects of the data: 
evolution in the color distribution of the massive galaxy 
population in (*}!]); the color evolution of the red galaxy 
population as a whole (37|); and the Zphot ^ 2 evolution 
in the absolute and relative numbers of massive red/blue 
galaxies (Sj8]). Our final results are in confiict with those 
from COMBO-17 in the same field; in Appendix [Bl we 
show that this is a product of calibration errors in the 
COMBO-17 data, rather than differences in our analyses. 
In SjH we present a series of sensitivity analyses in which 
we repeat our analysis a number of times, while varying 
individual aspects of our experimental design, and see- 
ing how these changes affect our results; this tests thus 
enable us to identify and quantify the most important 
sources of systematic uncertainty in our main results. 

Finally, in ijlOl we present a completely independent 
consistency check on our results: we measure the z < 2 
evolution of the relative number of bright, red galaxies 
based only on directly observed quantities — that is, 
without deriving redshifts or stellar masses for individual 
galaxies. A summary of our results and conclusions is 
given in Hill 

Throughout this work, magnitudes are expressed in 
the AB system; exceptions are explicitly marked. All 
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TABLE 1 
Summary of the data comprising 
THE MUSYC ECDFS catalog 



Band Aq Instrument Exp. Time Area FWHM 5a depth 









[mill.] 


rn'i 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


u 


3560 


WFI 


1315 


973 


I'.'l 


26.5 


U38 


3660 


WFI 


825 


942 


1"0 


26.0 


B 


4600 


WFI 


1157 


1011 


I'.'l 


26.9 


V 


5380 


WFI 


1743 


1020 


I'.'O 


26.6 


R 


6510 


WFI 


1461 


1016 


0'.'9 


26.3 


I 


8670 


WFI 


576 


976 


I'.'O 


24.8 


z' 


9070 


MOSAIC-II 


78 


997 


I'.'l 


24.0 


J 


12500 


ISPI 


Si 80 


882 


< 1'.'5 


23.3 


H 


16500 


Son 


60 


650 


< 0'.'8 


23.0 


K 


21300 


ISPI 


Si 60 


887 


< I'.'O 


22.5 



Note. — For each band, we give the filter identifier (Col. 1) 
and effective wavelength (Col. 2), as well as the detector name 
(Col. 3). For the NIR data, exposure times (Col. 4) are given per 
pointing; the effective seeing (Col. 6) is given for the pointing with 
the broadest PSF. The 5(t limiting depths (Col. 7) are as measured 
in 2'.'5 diameter apertures; the smallest apertures we use. The 
r efere nces for each set of imaging data are given in the main text 



masse s have been derived assuming a 'diet Salpeter' 
IMF ()Bell fc de Jongll2001h . whic h is defined to b e 0.15 
dex less massive than a standard ISalpeteil ()1955[ ) IMF. 
In terms of cosmology, we have assumed f^A = 0.70, 
fim = 0.30, and iJp = 70 /170 km s~^ Mpc~^, where 
hro = 1. 

2. DATA 

2.1. An Overview of the MUSYC ECDFS Dataset 

This work is based on a if-selected catalog of the 
ECDFS from the MUSYC wide NIR imaging pro- 
gramme; these data are described and presented in Paper 
I. We will refer hereafter to this dataset as 'the' MUSYC 
ECDFS catalog, although it should be distinguished from 
the optical (i?-|-F-|-i?)-selected catalog, and the narrow 
ba nd (5000 A)-selected photometric catalogs described 
bv lGronwall et all (l2007f). and th e spectroscopic catalog 
described bv lTreister et all (|2008D . 

The vital statistics of the imaging data that have gone 
into the MUSYC ECDFS catalog are summarized in Ta- 
ble □ Unlike the three other MUSYC fields, the ECDFS 
dataset was founded on existing, publicly available op- 
tical imaging: specifically, archival UUssBVRI WFI 
data, in cluding those tak en as part of the COMBO-17 
survey (iWolf et al.l 120031) and ESO's Deep Public Sur- 
vey (DPS: lArnouts et al.l |2001D , which h ave been re- 
reduced as part of the GaB oDS project (lErben et al.l 
120051: iHildebrandt et ani2006D . We also include H band 
imaging (P Barmby, priv. comm.) from SofI on the 3.6m 
NTT, covering ^ 80 % of the field, taken to c omplement 
the ESO DPS data, reduced and described bv lMov et al.l 

dioos). 

2.1.1. Original Data Reduction and Calibration 

These existing data have been supplemented with orig- 
inal z' band i maging from MOSAIC-II, reduced as per 
iGawiser et al.l ((2006), as well as J and K band imag- 
ing from ISPI; both instruments are mounted on the 



Blanco 4m telescope at CTIO. To cover the full 5 x 1 
ECDFS in the JK bands, we have constructed a mosaic 
of nine ~ 10 x 10 □' subfields (the size of the ISPI field 
of view ). The data reductio n for the JK imaging c losely 
follows iQuadri et al.l (|2007| ) and lBlanc erall (|2008) . and 
is described in detail in Paper I, where we present the 
MUSYC ECDFS catalog. 

In brief, to facilitate multiband photometry, each re- 
duced image has been shifted to a common astrometric 
reference frame (0'.'267 pix^^). The relative astrometry 
has been verified to O'.'IS (0.56 pix). To combat aperture 
effects (ie. similar apertures capturing different fractions 
of light, due to variable seeing across different images), 
we have PSF-matched our images to the one with the 
worst effective seeing. Among the K band pointings, the 
worst effective seeing is I'.'O FWHM; this sets our limits 
for detection and for total K band fiux measurements. 
Among the other bands, the worst seeing is 1'.'5 FWHM 
in the Eastern J subfield; this sets the limit for our multi- 
color photometry. After PSF matching, systematic errors 
due to aperture effects are estimated to be < 0.006 mag 
for the smallest apertures we use. 

We have tested the photometric calibration through 
comparison with the COMBO-17 catalog of the ECDFS 
(|Wolf et al.|[2003 l). and with the FIR EWORKS catalog 
of the GOODS-CDFS region (|Wuvts e t al. 2008). While 
there are significant differences between the COMBO- 
17 and MUSYC photometry, the comparison to FIRE- 
WORKS validates our photometry and photometric cal- 
ibration to < 0.02 mag in most cases, particularly for the 
redder bands. Further, we have tested the relative cal- 
ibration of all bands using the observed colors of stars; 
this test validates the photometric cross-calibration to 
< 0.05 mag. (See tj9.1l for a discussion of how sensitive 
our main results are to photometric calibration errors.) 

2.1.2. Photometry 

The photometry itse lf was done using SExtractor 
(|Bertin fc Arnoutdll996[) in dual image mode, using the 
I'.'O FWHM K mosaic as the detection image. Note that, 
as we were unable to find a combination of SExtrac- 
tor background estimation parameters (for the detection 
phase) that were fine enough to map real spatial varia- 
tions in the background level, but still coarse enough to 
avoid being influenced by the biggest, brightest sources, 
we were forced to perform our own background subtrac- 
tion for the NIR images. Total fluxes were measured from 
this I'.'O K image, using SExtractor's FLUX_AUTO. In 
Paper I, we that (in the photometry phase) SExtractor 
systematically overestimates the background flux level by 
~ 0.03 ma g; we have taken ste ps to correct for this effect. 
Following Labbe et al.l ()2003[ ). we also apply a minimal 
correction to account for missed ffux beyond the flnite 
AUTO aperture, treating each object as though it were 
a point source. We quantify the impact these two cor- 
rections have on our final results in i j9.2.1l and ij9.2.2l 

Multicolor spectral energy distributions (SEDs) were 
constructed for each object using the larger of SExtrac- 
tor's ISO aperture and a 2'.'5 diameter circular aperture, 
measured from the 1'.'5 FWHM U — K images; we then 
normalize each object's SED using the total K flux. This 
flexibility in aperture size is important to compromise 
between using apertures that are small enough to opti- 
mize S:N (the 2'.'5 diameter aperture is close to optimal 
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in terms of S:N for a point source in the 1'.'5 FWHM J 
band image), but also large enough to account for color 
gradients, which are important for the nearest, brightest 
objects. (See i )9.2.3l for a discussion of how our results 
vary using only fixed aperture photometry to construct 
SEDs.) 

Photometric errors (accounting for sky noise, imper- 
fect background subtraction, etc., as well as the pixel- 
pixel correlations introduced at various stages in the 
reduction process) were derived empirically by placing 
large numbers o f 'empty' apertures o n ea ch image (see 
als2 jLabbe et al.ll2003i iGawiser et aTll2006l : [Qu"adri et all 
|2007( ). For the J and K bands, this was done for each 
subfield individually. 

2.1.3. Completeness and Reliability 

We have assessed the completeness of the MUSYC cat- 
alog of the ECDFS in two ways (Paper I). First, we 
have tested our ability to recover synthetic, R^/^ pro- 
file sources of varying total flux and size introduced into 
empty regions of the data, using procedures identical to 
'live' detection. This analysis suggests that, at K = 22, 
the MUSYC catalog should be 100 % complete for point 
sources, dropping to 64 % for R^s — C.'G (2.25 pix, or 4.8 
kpc at z = 1) ellipticals. Secondly, we have compared our 
catalog to the much deeper FIREWORKS (|Wuvts et al.l 
l2008l ) catalog of the CDFS-GOODS region. In the region 
of overlap, for 21.8 < K < 22.0 bin, > 85 % of FIRE- 
WORKS detections are found in the MUSYC catalog; 
all MUSYC detections in this bin are confirmed in the 
FIREWORKS catalog. Taken together, these two anal- 
yses suggest that, at our limiting magnitude oi K — 22, 
the MUSYC catalog is primarily magnitude (c/. surface 
brightness) limited, > 85 % complete, and ^ 100 % reli- 
able. 

2.1.4. Sample Selection 

In constructing our main galaxy sample, we have iden- 
tified stars on the basis of the {B — z')-{z' — K) color- 
color diagram. This selection performs extremely well in 
comparison with both COMBO-17's SED classification, 
and with GEMS point sources (Paper I). We also make 
three further selections. First, to protect against false 
detections in regions with lower weights {e.g. the mosaic 
edges, and exposure 'holes' in the Eastern K subfield), 
we require the effective weight in K to be greater than 75 
%, corresponding to an effective exposure time of 45 min 
or more. Secondly, we have masked out, by hand, regions 
around bright stars where the SEDs of faint objects may 
be heavily contaminated; this problem is most severe in 
the z' band, where the PSF has broad wings. With these 
two selections, the effective area of the catalog becomes 
816 □'. Thirdly, to protect against extremely poorly con- 
strained redshift solutions, we will limit our analysis to 
those objects with S:N > 5 for the K SED point. (See 
ij9.3.1l for a discussion of how this selection impacts our 
results.) 

Of the 16910 objects in the MUSYC ECDFS catalog, 
these selections produce 10430 reliable i^T < 22 detec- 
tions, of which 9520 have reliable photometry in the ISPI, 
MOSAIC-II, and WFI band passes. Of these, 8790 cat- 
aloged objects have K S:N > 5; 950 of these objects 
are excluded as stars, leaving 7840 galaxies in our main 
sample. 



2.2. The z « Comparison Sample 

We will investigate the Zphot ^ 2 evolution in the mas- 
sive galaxy population by comparing the situation at 
z ^ to that for z ~ galaxies in the Sloan Digital 
Sky Survey fSDSS: r^riTetaniMol) : specifically, we use 
the 'low-z' sample from the New York University (NYU) 
Value Ad ded Galaxy Ca t alog (V AGC) of the SDSS pre- 
sented bv lBlanton et all ()200 5b\ The (Data Release 4) 
low-z catalog contains ugriz photometry for 49968 galax- 
ies with IQ < D < 150 Mpc (zspcc < 0.05), covering 
an effective area of 6670 □"; 2513 of these galaxies have 
Af, > 10^^ M0. Our analysis of these data closely follows 
that of the z ;2> sample, and is described separately in 
Appendix [X] 

3. PHOTOMETRIC REDSHIFTS 
AND RESTFRAME PROPERTIES 

3.1. Photometric Redshifts 

The technical crux on which any photometric look- 
back study rests is the determination of redshifts from 
broadband SEDs. We h ave computed our ph otomet- 
ric redshifts using EAZY H rammer "eTal] 120081 ). a new, 
fully-featured, user-friendly, and publicly-available pho- 
tometric redshift code. By default, the Zphot calculation 
is based on all ten bands, although we do require that the 
effective weight in any given band is greater than 0.6; in 
practice, this requirement only affects the H band, where 
we do not have full coverage of the field. The character- 
istic filter response curves we use account for both at- 
mospheric extinction and CCD response efhciancy as a 
function of wavelength. 

For our fiducial or default analysis, we simply adopt 
the recommended default settings for EAZY: viz., we 
adopt the EAZY default redshift and wavelength grids, 
template library, template combination method, tem- 
plate error function , K luminosity prior, etc. (See 
iBrammer et al.l [20081 for a complete description.) Note 
that, by default, EAZY assigns each object a redshift by 
marginalizing over the full probability distribution rather 
than, say, through minimization. For this work, a key 
feature of EAZY is the control it offers over how the SED 
fitting is done: the user is able to specify whether and 
how the basic template spectra are combined, whether or 
not to include luminosity prior and/or a template error 
function, and how the output redshift is chosen. We will 
make use of EAZY's versatility in ^9.4\ to explore how 
these particular choices affect our results. 

One of the unique aspects of the ECDFS is the high 
number of publicly available spectroscopic redshift de- 
terminations, which can be used to validate and/or cali- 
brate our photometric redshifts. In Paper I, we describe a 
compilation of spectroscopic redshifts for 2914 unique ob- 
jects in our catalog, including 'robust' redshifts for 1656 
galaxies in our main K < 22 sample. These redshifts 
come from some of the many literature sources available 
in the ECDFS, including those large surveys referred to 
in SJTJ as well a s the X-ray selected spec troscopic red- 
shift c atalogs of ISzokolv et al.l (I2004D and Treister et all 
(f20Q8|), a new survev bv iS Koposov et al.l (jin prep.l ). and 
a number of smaller projects. K20 is particularly use- 
ful in this regard, given its exceptionally high spectro- 
scopic completeness, albeit over a very small area: 92 % 
of A'(^'=s'') < 20 sources over 52 □'. 
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Fig. 1. — Validating the MUSYC catalog photometric redshift 
determinations — Main panel: the ^spcc^-Zpijot diagram for the 1656 
galaxies from our main K < 22 sample, using a compendium of 
'robust' spectroscopic redshifts from the literature (see Paper I). 
Inset: the distribution of Az = Zphot — ^spcc for the same set of 
galaxies; the curve shows a Gaussian fit to this distribution, with 
parameters as given. Lower panel: the redshift distributions of 
the main and spec-z samples; our photometric redshifts appear 
to mildly underestimate the redshift of the three overdensities at 
0.5 < Zspec < 0.8. Quantitatively, we find the median and NMAD 
of Az/{1 + Zspec) for the full spec-z sample to be -0.029 and 0.036, 
respectively (See also Figure^; for Zspec > 1, we find these num- 
bers to be -0.023 and 0.060; for the K20 sample, which is 92 % 
complete for /ft^B) ^ 2I.8, these numbers are -0.028 and 0.033. 



The main panel of Figure [T] shows our Zspcc-^phot plot. 
We prefer to quantify the photometric redshift quality 
in terms of the normalized median absolute deviation 
(NMAD^) in Az/(1 -I- ^spoc), which we will abbreviate as 
az', for this comparison sample, az — 0.035. Further, the 
outher fraction is acceptably small: 5.9 %. Comparing 
only to the 241 redshifts from K2 0, we find Uz = 0.033; 
for the iVan der Wei et all (|2005 f) sample of 28 z ~ 1 
early type galaxies the figure is 0.022. For 1 < Zgpcc < 2, 
we find az = 0.059. For the 20 % (269/1297) of 0.2 < 
Zphot < 1.8 galaxies from our mass limited sample defined 
in i|4] that have spectroscopic redshifts, we find az = 
0.043. 

Ba sed on their catalog of the GOODS ACS and ISAAC 
data, iGrazian et al.l ()2006( ) have achieved a photometric 
redshift ac curacy of (Az/ ( 1 + Zs pon)) of 0.045. For com 



parison to IGrazian et alj ()2006( ). the inset panel shows 

^ The NMAD is defined as 1.48 X median[|a; — median(a;)|]; the 
normalization factor of 1.48 ensures that the NMAD of a Gaussian 
distribution is equal to its standard deviation. 



the distribution of Az = (zphot — Zspcc) for < Zspoc < 2; 
although offset by -0.046, the best fit Gaussian to the 
dist ribution has a width of 0.046, as opposed to 0.06 
for IGrazian et al.l ()2006D . For an identical sample of 
938 galaxies with Z spccSi we find az — 0.043 for the 
IGrazian et al.l ()2006f ) 2:photS and az = 0.035 for ours. 
In other words, our photometric redshift determinations 
are at least as good as the best published for iiT-selected 
samples at high redshifts. We also note in passing that 
our Zphot ^ 1 photometric redshifts agree very well with 
those from COMBO-17 (Wolf e t al.. ,20031 : a detailed 
comparison to both these catalogs is presented in Paper 
I. 

The lower panel of Figure [T] shows the redshift dis- 
tributions for both our main galaxy sample (based on 
•Zphot), and the spectroscopic comparison sample (based 
on Zspec)- Note the presence of three prominent red- 
shift spikes at 0.5 < Zspec < 0.8 (see also IVanzella et ahl 
l2008h : it appears that our redshift determinations may 
slightly underestimate the redshifts of these structures. 
The structures at z ~ 1.0, 1.1, 1.2, 1.3, and 1.4 
(jVanzella et al.l [20081 ) are also visible in the Zgpec distri- 
bution, but are 'washed out' in the Zphot distribution. 

3.2. Restframe Photometry and Stellar Masses 

The many degeneracies between SED shape and the 
intrinsic properties of the underlying stellar population, 
which are actually a help when deriving photometric red- 
shifts, make the estimation of such properties from SED 
fitting highly problematic. Systematic uncertainties as- 
sociated with parameterisations of the ass umed star for- 
matio n history are at the level of 0.1 dex (jPozzetti et al.l 
l2007f ). while uncertainties in the stellar population mod- 
els themselves are generally accepted to be < 0.3 dex; 
this is comparable to the uncertainty associated with the 
choice of stellar IMF. For these reasons, we have opted 
for considerably simpler means of deriving restframe pa- 
rameters. 

Once the redshift is determined, we have interpolated 
restframe fluxes from the observed SED using a new util- 
ity dubbed InterRest, which is a slightly more sophisti- 
ca ted version of the alg orithm described in Appendix C 
of lRud nick et al.l (|2003f l. and is described in detail in Pa- 
per I. InterRest is designed to dovetail with EAZY, and 
is also freely available.^ We estimate the systematic er- 
rors in our interpolated fluxes (c/. colors) to be less than 
2 % (Paper I). 

We then use this interpolated restframe photometry 
to es timate galaxies' stellar masses using a prescription 
from Bell &^de_Joiig| (|2001f) . which is a simple linear re- 
lation between restframe {B — V) color and stellar mass- 
to- light ratio: M^,/Lv- 



logio M^/Li 



.734 + 1.404 X (B-y + 0.084) , (1) 



assuming My.o ~ 4.82. (Here, the factor of 0.084 is 
to transform fr o m the Vega magnitude system used by 
iBeh fc de Jond (|2001f) to the AB system used in this 
work.) This prescription assumes a 'diet Salpeter' IMF, 
which is defined t o be 0.15 dex less massive than the 
standard lSalpeteil (|1955D IMF, and is approxim ately 0.04 
dex more massive than that of iKroupal (j2001f ). Further, 

Code and documentation can be found at: 
|http: //www. strw.leidenuniv.nl/~ent/Inter Rest/ 1 



6 



E N Taylor et al. 



to prevent the most egregious overestiraates of stellar 
masses, we limit M^jLy < 10 (see also Figure 4 of Borch 
et al. 2006). Although this Hmit affects just 1.2 % of our 
main sample, we found it to be important for getting 
the high- mass end of the z w mass function right (Ap- 
pendix [A| . 

It is not immediately obvious whether using these 
color-derived M/Ls is significantly better or worse than, 
say, from using stellar population synthesis to fit the 
whole observed SED. The prescription we use has been 
derived from SED-fit M/Ls; the scatter around this re- 
lation is on the order of 0.1 dex. By comparison, the 
precision of SED-fit M/L determinations is limited to 
0.2 dex by degeneracies between, e.g., age, metallic- 
ity, and dust obscura tion (see, e.g.. iPozzetti et aLll2007l : 
iConsehce et~aIll2007D . Thus, the increase in the random 
error in M^, due to the use of color-derived rather than 
SED-fit stellar masses is ~ 10 %. 

In addition to being both simpler and more transpar- 
ent, however, the use of color-derived M/Ls has the ma- 
jor advantage of using the same restframe information 
for all galaxies, irrespective of redshift. This is especially 
important when it comes to the comparison between the 
high-z and low-z samples, where the available photome- 
try samples quite different regions of the restframe spec- 
trum. Since Equation[l]has ultimately been derived from 
SED-fit mass estimates, however, our color-derived mass 
estimates are still subject to the same systematic uncer- 
tainties. To the extent that such systematic effects are 
independent of color, redshift, etc., they can be accom- 
modated within our results by simply scaling our limiting 
mass. On the other hand, if there is significant evolution 
in the color-M/L relation with redshift, then there is the 
risk that the use of color-derived M/Ls may introduce se- 
rious systematic errors with redshift. We investigate this 
issue further in §9.5) in which we also demonstrate that 
our results are essentially unchanged if we use conven- 
tional SED-fitting techniques to derive M/Ls. 

3.3. The Propagation of Redshift Errors 

A primary concern in this paper is the importance of 
systematic errors. To address this concern in the context 
of our photometric redshift determinations, we show in 
the top panels of Figure O the 2;spec--Zphot agreement as 
a function of redshift, S:N in the K 'color' aperture, and 
restframe color. In each panel, the points with error bars 
show the median and 15/85 percentiles in discrete bins. 

Looking first at our photometric redshifts: the first 
panel of Figure [H shows the photometric redshift error, 
Az/ (l-l-ZspGc), as a function of Zgpec- We see that there is 
a systematic offset between Zphot and ^spec for Zgpec ^ 1, 
such that our ZphotS tend to be slightly too low (see also 
Figure [1]); for z > 1, this effect appears to be less. At 
least for K S:N > 10, random errors in the photometric 
redshifts do not appear to be a strong function of S:N. 

There is a clear systematic effect as a function of rest- 
frame color. For galaxies redder than {u — r) « 2 (ap- 
proximately the lower limit for z « red sequence galax- 
ies), the agreement between Zphot and Zgpoc is very good. 
For galaxies with (u — r) < 2, however, it seems that 
we systematically underestimate the true redshift by ap- 
proximately Az < 0.02(1 -I- Zgpcc). It is plausible that 
this is in fact the driver of the weak apparent systematic 
effect with redshift, coupled with there being a greater 



proportion of blue galaxies in the spectroscopic redshift 
sample at lower redshifts. 

How do these errors in redshift estimation play out 
in the derivation of restframe properties? We address 
this question with reference to the lower panels of Fig- 
ure [2l which illustrate how redshift uncertainties affect 
our derivation of three basic restframe quantities {top to 
bottom): absolute luminosity, restframe color, and stellar 
mass. In each panel, we plot the difference between the 
values derived adopting the spectroscopic or photometric 
redshift, as a function of {left to right) redshift, K band 
S:N, and restframe color, as well as photometric redshift 
error. 

Quantitatively, for our Zgpcc comparison sample, the 
random photometric redshift error of Az/(1 + Zgpcc) = 
0.035 translates into a 0.360 mag error in absolute mag- 
nitude, 0.134 mag error in restframe color, and a 0.107 
dex error in stellar mass. (By contrast, the typical uncer- 
tainty in for a X S:N = 10 galaxy is AK = 0.12-0.16 mag 
« 0.05-0.07 dex.) Just as for the redshifts themselves, 
the clearest systematic effects in the derived quantities is 
with restframe color: there appear to be mild systemat- 
ics with redshift for z < 1, and no clear trend with S:N, 
at least for S:N > 10. 

It is straightforward to understand why redshift errors 
play a larger role in the derivation of magnitudes rather 
than colors. When calculating magnitudes, the primary 
importance of the redshift is a distance indicator. For 
the Zgpoc sample shown in Figure [21 the random scatter 
in AA'Ir due to distance errors alone (calculated by taking 
the difference in the distance modulus implied by Zphot 
versus that by Zgpcc) is 0.28 mag; ie. ^ 75% oi the scatter 
seen in Figure [21 On the other hand, colors are distance 
independent, and this clement of uncertainty is cancelled 
out. 

What is surprising is the relative insensitivity of our 
stellar mass estimates to redshift errors. Focusing on 
the right panels of Figure [2l it can be seen that where 
the photometric redshift underestimates the true red- 
shift/distance, we will infer both too faint an absolute 
luminosity and too red a restframe color. When it comes 
to computing a stellar mass, however, these two effects 
operate in opposite directions: although the luminosity 
is underestimated, the too-red color leads to an overes- 
timate of the stellar mass-to-light ratio. The two effects 
thus partially cancel one another, leaving stellar mass 
estimates relatively robust to redshift errors. 

In a photometric redshift survey, the measurement un- 
certainty due to random photometric redshift errors is 
considerably less for stellar masses than it is for absolute 
magnitudes. This conclusion remains unchanged using 
SED-fit stellar masses, rather than our favored color- 
derived ones. Conversely, we can say that (random) 
photometric redshift errors are not a dominant source of 
uncertainty in our stellar mass estimates. Indeed, as we 
have already noted, the size of these errors is comparable 
to the uncertainties in our total flux measurements. 

4. CONSTRUCTING A 
STELLAR MASS SELECTED SAMPLE 

For moderate to high redshifts, NIR selection has the 
key advantage of probing the restframe optical light, 
which is a reasonable tracer of stellar mass. In this 
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Fig. 2. — Photometric redshift errors, and their effect on other derived quantities — In each panel, the abscissa shows the difference in 
a derived quantity, derived assuming the spectroscopic or photometric redshift; in all cases, 'A' should be understood as the difference 
between the Zphot^ minus ^spcc— derived values. We show: (top to bottom, in rows) redshift, absolute magnitude, restframe color, and 
stellar mass, as a function of (left to right, in columns) redshift, observed signal— to-noise, restframe color, and photometric redshift error. 
The black points are for the spectroscopic sample shown in Figure [T] the red points show the median offset in bins; the error bars reflect 
the 15/85 percentiles. In each panel, points that fall outside of the plotted range are shown as small grey plusses. The separate panels at 
right show the distribution of As for all galaxies in the Zspcc sample; the median and scatter (NMAD) in the difference between Zphot^ a^nd 
Zspcc-dcrivcd quantities are as shown. In all cases, the clearest systematic effects are as a function of redshift, with some systematic errors 
for the reddest and bluest galaxies. The way that redshift errors propagate mean that the uncertainties due to redshift errors are much 
smaller for stellar masses than they are for either absolute magnitudes or restframe colors. 



section, we empirically relate our observed flux detec- 
tion/selection limit to an approximate completeness limit 
in terms of stellar mass and redshift. 

To this end, we have taken galaxies with K fluxes im- 
mediately below our detection limit from three signifi- 
cantly deeper f T-selected catalogs; tizz. the MUSYC deep 
NIR catalogs ("Quadri et alj [20071 ). the FIREWORKS 
catalog dW uvts ct al. 2008 ), and the FIRES ca talogs 
(|Labbe et al.l i2003: Fo rster-Schreiber et all I2006D . By 
taking objects from these catalogs that lie immediately 
below our detection threshold, and scaling their fluxes 
(and so stellar masses) to match our K = 22 limit, it is 
then possible to empirically determine the stellar mass- 
redshift relation for K fa 22 galaxies. The upper envelope 
of points in (M„, Zphot) space thus represents the most 
massive galaxies at our observed flux limit, and so di- 
rectly provides a redshift-dependent mass completeness 
limit. 

This is illustrated in the left panel of Figure [31 In 



this panel, the large, open, colored symbols represent 
22.0 < K < 22.5 objects from the deeper catalogs, scaled 
up in flux to = 22; viz. the MUSYC deep NIR cat- 
alogs {blue crosses), the FIREWORKS catalog {yellow 
circles), and the FIRES catalogs {red squares). Again, 
these points represent objects immediately at our detec- 
tion limit; the upper envelope of these points therefore 
represents the most massive galaxies that might escape 
detection/selection in our analysis. This suggests that 
for Af, > 10^^ Mq, we are approximately complete for 
2phot < 1-8. 

It is possible to do the same thing using the faintest 
detections from our own catalog, scaled down in flux to 
our selection limit. Specifically, we have taken galax- 
ies with 21.5 < K < 22.0 and scaled their fiuxes (and 
masses) down to i^T = 22. For this test, we also restrict 
our attention to galaxies with well constrained redshifts, 
by requiring that the EAZY 'odds' parameter be greater 
than 0.95. 
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Fig. 3. — Empirically determining our mass completeness limit as a function of redshift — Left panel — The black points show stellar 
masses for MUSYC ECDFS galaxies with 21.5 < K < 22.0, scaled down in flux to match our K = 22 detection limit, and plotted as 
a function of photometric redshift. The other symbols show stellar masses for 22.5 > K > 22.0 galaxies, scaled up in flux to K = 22; 
these galaxies are drawn from the MUSYC deep fields (blue crosses), the FIREWORKS catalog (yellow circles), and the FIRES catalogs 
(red squares). Each sample has been analyzed in exactly the same manner. The upper envelope of these points effectively defines, as a 
function of redshift, the limiting stellar mass corresponding to our observed K flux limit. For Af* > lO^'^ Mq, we are nearly complete 
(> 90%) to 

^phot — 1-8. Right panel — the color— stellar mass diagram for K ^ 22 galaxies at 1.7 < 2phot < 1-9 — the large squares show the 
median values for the MUSYC ECDFS points, binned by color; all other symbols are as in the left panel. Here, the right envelope of the 
colored points defines our mass completeness limit at Zphot ~ 1.8 as a function of color. For comparison with Figures |4] and (5] the hatched 
region shows estimated completeness limits based on synthetic SSP spectra. While we may well miss galaxies considerably redder than the 
predicted red sequence (see JTJl, at Zphot ~ l-Sj this empirical argument suggests that we are approximately complete (> 85 %) for galaxies 
with Af, > 10" M0 and (u - r) < 2. 



These points are shown as the closed circles in the left 
panel of Figure[21 While the results of this 'internal' test 
are broadly consistent with the previous 'external' one, 
they do suggest slightly higher incompleteness. Of the 
21.5 < K < 22.0 sources with 1.6 < Zphot < 1-8, 23 % 
(20/87) would have > 10" M© when scaled down to 
K = 22, indicating that our completeness for K = 22, 
= 10" M0 galaxies is ^ 75 % for 1.6 < Zphot < 

1.8. However, the 21.5 < K < 22.0 subsample shown 
here represents only 30 % of our full K < 22 sample in 
this mass and redshift range, suggesting that the overall 
completeness is more like > 90 %. 

As a second and complimentary check on this conclu- 
sion, the right panel of Figure [3] shows the color-stellar 
mass diagram for a narrow redshift slice at 1.7 < Zphot < 

1.9. Here, the large squares show the median values from 
the MUSYC ECDFS points, binned by color; the other 
symbols are the same as in the other panel of this Fig- 
ure. As before, the points in this panel represent ob- 
jects at our detection limit; the right envelope of these 
points thus describes our mass completeness limit for 
•Zphot ~ 1-8, this time as a function of restframe color. 

Again, the 'internal' and 'external' analyses broadly 
agree. For blue galaxies, both tests suggest that MUSYC 
should be very nearly complete for M^, > 10^^ Mq and 
Zphot < 1-8. For {u — r) > 1.5 galaxies, however, the 
down-scaled MUSYC points again suggest slightly lower 
completeness than those scaled up from deeper catalogs: 
45 % (13/29) of these galaxies would fall foul of one of 
our selection criteria if their masses were scaled down to 



= 10^"'^ Mq. Using an argument analogous to that 
above, this suggests that our completeness fraction for 
galaxies with M* > 10-^^ Mq, Zphot = 1-8, and (u — r) > 
1.5 is at least ~ 85 %. 

We therefore adopt Af* > 10" Mq and Zphot < 1-8 
as our approximate completeness limits, corresponding 
to our K < 22 detection/selection limit. In ii9.3.2i we 
apply simple completeness corrections to determine the 
extent to which our results may be affected by incom- 
pleteness. As a final caveat, however, there remains the 
concern of additional incompleteness due to our S:N 
> 5 criterion, which we will address in ij9.3.1l 

5. THE COLOR-MAGNITUDE AND 
COLOR-STELLAR MASS DIAGRAMS FOR Zphot < 2 

In this section, we present our basic observational re- 
sults: the color-magnitude and color-stellar mass dia- 
grams for Zphot < 2. 

5.1. The Color-Magnitude Diagram for Zphot ^ 2 

Figure m shows the color-magnitude diagram (CMD), 
plotted in terms of {u—r) color and absolute r magnitude, 
Mr, for Zphot < 2. 

The first panel of Figure[5]is for z « galaxies from the 
'low-z' comparison sample; these data and our analysis 
of them are described in Appendix [Xl The basic features 
of the CMD — the red sequence, blue cloud, and green 
valley — are all easily discernible. The dotted line shows 
our characterization of the z CMR for red galaxies 
(Equation [XT]), also discussed in Appendix [XI 
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Fig. 4. — The color— magnitude diagram (CMD) for galaxies with Zphot ^ 2 — The first panel shows a random selection from the NYU 
VAGC's 'low-z' sample, based on DR4 of the SDSS (Blant on et al.ir2005bl V discussed in Appendix lAl the other panels show the MUSYC 
ECDFS data, discussed in the main text. Except where marked, bins are of equal comoving volume; for the low-z sample, we have plotted a 
random sub-sample to yield the same effective volume: the density of points is thus directly related to bivariate comoving number density. 
The shaded area shows the approximate K = 22 detection/selection limits, based on synthetic spectra for an SSP; the error bars show 
representative errors for a M« 10^^ Mq, (m — r) « 2.0 galaxy at the mean redshift of each bin. The dotted line in each panel shows our 
fit to the OMR for bright, red sequence galaxies at z ^ 0, derived in Appendix |A] In this work, we prefer to use stellar mass, rather than 
absolute magnitude, as a basic parameter — accordingly, we focus our attention on the CM.D, presented in Figure[5] 



The other eight panels show the 0.2 < Zphot < 1-8 
MUSYC ECDFS data. Except for the two highest red- 
shift bins, which are twice as large, the Zphot ^ bins 
have been chosen to have equal comoving volume.^ For 
the z w bin, we plot only a random sub-sample of the 
low-z catalog, chosen to effectively match the volume of 
the higher redshift bins. Thus, the density of points in 
the color-magnitude plane is directly related to changes 
in the bivariate comoving number density. 

In the bottom-right corner of each panel, we show rep- 

^ The exact redshift limits we have used are ^phot = 0.200, 0.609, 
0.825, 0.987, 1.127, 1.254, 1.373, 1.486, 1.595, 1.700, 1.804, 1.906, 
and 2.000; elsewhere we will round these values as convenient. 



resentative error bars for a ss 10"'^"'^ Mq galaxy with 
(u — r) « 2.0, near the mean redshift of each bin. In 
order to derive these errors, we have created 100 Monte 
Carlo realizations of our catalog, in which we have per- 
turbed the catalog photometry according to the photo- 
metric errors, and repeated our analysis for each: the 
error bars show the scatter in the values so derived. The 
shaded grey regions show approximately how our K < 22 
completeness limit projects onto color-magnitude space 
through each redshift bin, derived using synthetic single 
stellar population (SSP) spectra. 

Examining this diagram it is clear that, in the most 
general terms possible, bright/massive galaxies were con- 
siderably bluer in the past. At a fixed magnitude, the 
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Fig. 5. — The color— stellar mass diagram (CM,D) for galaxies with Zphot ^2 — As in Figure |4] the z bin is based on the low-z 
sample of SDSS galaxies, discussed in Appendix |X] the Zphot ^ points based on the MUSYC ECDFS data, described in the main text. 
The hatched area shows approximate selection limits, based on synthetic spectra for an SSP; our empirical completeness limit is marked 
in the last bin. The error bars show representative errors for a M, 10" Mq, {u-r) ^ 2.0 galaxy near the bin's mean redshift, based on 
100 Monte Carlo realizations of the catalog data, including photometric redshift errors. Within each panel, the dotted line shows our fit to 
the CM»R for bright red sequence galaxies at z 0, derived as per Appendix|A] for the Zp^ot 2> bins, the solid lines show our fit to the 
color evolution of the massive red galaxy population, derived in JT] the dashed line shows our red galaxy selection criterion, introduced in 
3H] We analyse the key features of this diagram further in Figures |H] |5] O and IIUI 



entire z ~ 1 galaxy population is a few tenths of a magni- 
tude bluer than at z « 0. At the same time, particularly 
for z > 1, there is a growing population of galaxies with 
Mr < —22 and {u — r) < 2 that has no local analogue. 
While there are some indications of a red sequence within 
the z 3> data, particularly for Zp^ot 1> it is certainly 
not so easily distinguishable as locally. 

5.2. The Color-Stellar Mass Diagram for Zphot ^ 2 

There are a number of advantages to using stellar mass 
as a basic parameter, rather than absolute magnitude. 
Principal among these is the fact that stellar mass is 
more directly linked to a galaxy's growth and/or assem- 
bly: while a galaxy's brightness will wax and wane with 



successive star formation episodes, a galaxy's evolution 
in stellar mass is more nearly monotonic. On the other 
hand, it must be rememebered that the necessary as- 
sumptions in the derivation of stellar mass estimates pro- 
duce greater systematic uncertainties than for absolute 
luminosities. As bursts of star formation and other as- 
pects of differing star formation histories among galaxies 
are likely greater at higher redshifts, however, we will 
focus on the color-stellar mass diagram in this and fol- 
lowing sections. 

Figure [5] shows the color-stellar mass diagram (CM*D) 
for z < 2. As in Figure [H the first panel shows a random 
sub-sample of the low-z sample; the other panels show 
the MUSYC ECDFS data. The dotted line in each panel 
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Fig. 6. — Color distributions for M« > 10^^ Mq galaxies at 2 < 2 — In each panel, the histograms show the color distribution for 
Mt > 10^^ M0 galaxies, after subtracting out the slope of the local CMtR, normalized at M* = lO'^^ M0. The shaded distributions show 
the data themselves. The smooth curves show double Gaussian fits to the observed distributions. For 2phot ^ 1-2, the Gaussian fits to the 
observed color distributions are not robust; our inability to reliably distinguish separate red and blue populations for Zphot ^ 1-2 is likely 
due to insufficient S:N in our NIR data (see Figure[7]l. In each panel, the vertical dotted line shows the location of the z x red sequence; 
the vertical solid lines show our fit to the observed color evolution of the red sequence, derived in [[T] the dashed lines show our red galaxy 
selection criterion, introduced in flB] 



shows the z « color-stellar mass relation (CM*R), 
which we have derived in Appendix [Xj given in Equa- 
tion EH 

Each of the basic features of the CMD are also seen 
in the CM*D. We see an increasing number of galaxies 
with Zphot > 1 and with > 10" Mq and (u-r) <2 
which have no analogues in the local universe. For a SSP, 
the colors of these galaxies would suggest ages of < 1 
Gyr: these massive galaxies appear to be in the throes 
of their final star formation episodes. For Zphot ?J 1-2, 
these galaxies may even dominate the massive galaxy 
population. Some evidence for a distinct red sequence is 
visible in the CM*D for Zphot ^ 1, but not much beyond. 

The next three sections are devoted to more quanti- 
tative discussion of each of the following three specific 



observations: 

1. We see evidence for a red galaxy sequence for 
.Zphot ^ 1-2; beyond this redshift, whether due to 
physical evolution or to observational errors, it be- 
comes impossible to unambiguously identify a dis- 
tinct red sequence on the basis of the present data 

2. At a fixed mass, a red sequence galaxy at z ~ 1 is 
a few tenths of a magnitude bluer than its z ~ 
counterpart 

3. At higher redshifts, there appear to be fewer mas- 
sive galaxies on the red sequence. Further, it ap- 
pears that the proportion of blue cloud galaxies 
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among the most massive galaxies increases; con- 
versely, the red fraction is lower at higher redshifts. 

m- 

6. THE COLOR DISTRIBUTION 
OF MASSIVE GALAXIES FOR Zphot < 2 

In an attempt to quantitatively separate the mas- 
sive galaxy population into distinct red and blue sub- 
populations, Figure [6] plots the color distribution of the 
Af, > 10^^ Mq galaxy population, after subtracting out 
the slope of the z « CM,R, and using the same red- 
shift bins as in Figures |4] and \5\ The grey histograms 
in each panel show the data themselves. Note that for 
the z « panel of this plot, we have used the full low-z 
sample. Also, recall that for this mass regime, we are 
approximately complete (volume limited) to Zphot = 1-8. 

6.1. The Massive Galaxy Red Sequence at Zphot ^ 1-2 

Locally, red sequence galaxies totally dominate the 
massive galaxy population: what bimodality exists be- 
tween the red and blue populations is weak. (As the 
name suggests, 'bimodality' implies two distinct local 
maxima in the distribution.) This is a reflection of the 
apparent 'trans ition mass' betwe e n red and blue galax- 
ies observed bv lKauffmann et ahl ([2Q03); the bimodality 
is stronger in a luminosity limited sample, including a 
greater proportion of bright but less massive blue cloud 
galaxies (see also Figure \Al\ Appendix |A| . At slightly 
higher redshifts, where some progenitors of z w red se- 
quence galaxies are still forming stars in the blue cloud, 
we may then expect the bimodality to actually become 
stronger, before weakening again as the fraction of those 
galaxies already on the red sequence becomes small at 
moderate-to-high redshifts. In general, however, the 
color distributions shown in Figure [5] are not clearly bi- 
modal. 

With this in mind, as a simple means of separating 
red from blue galaxies, we have fit the observed distribu- 
tions in each redshift bin with double Gaussian functions. 
These fits are shown by the smooth curves in each panel 
of Figure [6l For the most part, these fits provide a rea- 
sonable description of the 0.2 < Zphot ^ 1-2 data (see 
also [Borch et al.ll2006l) . 

6.2. A Massive Galaxy Red Sequence at Zphot ^ 1.3? 

In contrast to lower redshifts, for Zphot ^ 1-3, we are no 
longer able to reliably fit the color distributions in this 
manner: on the basis of the sensitivity tests presented 
in ij9l neither the red/blue separation nor the fits to the 
distributions of these populations is robust. Looking at 
the typical errors shown in each panel of Figure [5l the 
measurement errors on the (u — r) colors of red galaxies 
rises sharply from A(-u — r) < 0.1 mag for Zphot ^ 1-2 
to 0.2 mag for Zphot ~ 1-3, 0.3 mag for Zphot ~ 1-5, and 
0.4 mag for Zphot ~ 1-7. This would suggest that our 
inability to reliably distinguish both a red and a blue 
population at these higher redshifts may very well be due 
to the large errors in restframc colors for higher redshift 
galaxies. 

In order to help interpret our Zphot ^ 1 results, we have 
tested our ability to recover a single-color galaxy popu- 
lation, given our observational and analytical errors. To 
this end, we have generated a mock photometric catalog 



containing only red sequence galaxies: beginning with 
our main galaxy sample, we have replaced each galaxy's 
photometry using synthetic spectra for a SSP formed 
over a short period beginning at Zform — 5, assuming 
the catalog values of Zphot and M*, and then adding 
typical MUSYC ECDFS photometric errors. We have 
then re-analyzed this catalog using the same methods 
and procedures as for the main analysis, including re- 
computing photometric redshifts, restframe colors, and 
stellar masses. 

The results of these tests are shown in Figure [3 which 
plots the observed width of the color distribution for this 
intrinsically single color population, as a function of pho- 
tometric redshifts, assuming photometric errors typical 
for the MUSYC ECDFS (squares). The measurement er- 
rors on the (u — r) colors of red galaxies rise sharply from 
0.05 — 0.07 mag for Zphot ^ 1 to 0.10 mag for Zphot ~ 1, 
and then continue to increase for higher redshifts. In or- 
der to demonstrate that this is a product of photometric 
errors per se and not redshift errors, we have also re- 
peated this analysis holding the redshifts of each object 
fixed; the results of this test are shown as the dotted line. 
Even with spectroscopic redshifts, the depth of our NIR 
data would seem to preclude the detection of a distinct 
red sequence at z > 1.3. 

What then would be required in order to confirm the 
non/existence of a red sequence at Zphot ^ 1-5? We 
have also constructed mock galaxy catalogs with S:N 
that is 2.5 and 10 times greater than typical values for 
the MUSYC ECDFS catalog; the results of these tests 
are shown as the triangles and circles, respectively. Even 
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Fig. 7. — The effects of photometric errors on our ability to 
recover a (red) single— color galaxy population — These results are 
based on mock galaxy catalogs containing only passively evolving 
galaxies, taking redshifts and stellar masses of galaxies from the 
main sample, which we have analyzed in the same manner as the 
actual data: we are thus testing our ability to recover a single-color 
galaxy population as a function of redshift. The squares, triangles, 
and circles show the width of the recovered color distributions as- 
suming typical errors for the MUSYC ECDFS data divided by 1, 
2.5, and 10, respectively; the dotted line shows the results assuming 
MUSYC ECDFS photometric errors, but no redshift errors. Even 
with spectroscopic redshifts, the depth of the MUSYC ECDFS NIR 
data precludes the detection of a Zphot ^1-4 red sequence. 
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Fig. 8. — The color evolution of massive galaxies for z <^ 2 — 
Points show the fit centers of the color distributions for the red and 
blue galaxy subpopulations (see Figure |6ll. Only the Zpijot ^1-2 
points {circles with error bars) were used when fitting for the color 
evolution of the red galaxy sequence {solid line); the errors were 
derived from bootstrap resampling. The large square at 2 = 2 
shows the approximate equivalent {u—r) color of the 3.3cr detection 
of a { U — B) red sequence among 2^2 galaxies from Kriek et al. 
1I2OO8I) . based on NIR spectra. We see rather a rather smooth 
reddening of the red sequence from 2 ~ 1.2 to the present day, 
which is well described by the linear fit given. 



Fig. 9. — Comparing the color evolution of massive galaxies for 
2 < 2 to passively evolving stellar population models — The data 
in this panel are identical to Figure |8] The overlaid curves show 
simple expectations for the passive evolution of a SSP, formed in 
a short burst (e-folding time of 100 Myr) beginning at 2foi.m = 2, 
3, 5, and 10 {bottom to top). While the Zf^Tm ^ 3 curves pro- 
vide a good qualitative description of the observed evolution, they 
have substantially lower metallicity than would be expected from 
the local mass-metallicity relation. Moreover, these very simple 
models make no attempt to account for progenitor bias, and other 
important effects. 



pushing a full magnitude deeper, it would be difficult to 
identify a red sequence at z 1.5, assuming that obser- 
vational errors of A(u — r) < 0.1 mag would be required 
to robustly identify a red sequence. In order to probe 
z > 1.5, an order of magnitude improvement is required. 
This would suggest that the detection of a red sequence 
at z > 1.5 would require a J band (5(7 point source) 
limit of ~ 25.8, roughly the final target depth for the 
Ultra D eep component of the UKIRT Infrared Deep Sky 
Survey ([Lawrence et al.l[2007f ). 

It is clear from this analysis that we cannot confirm or 
exclude the existence of a red galaxy sequence at z > 1.3 
on the basis on the present data. This is in good ac- 
cord with t he recent detections of a z < 1.5 red galaxy 
sequence by ICirasuolo et al.l (|2007( ) and Franzetti et al. 
(2007 — see also Kriek et al. 2008). Moreover, we note 
in passing that if we were to subtract away the broaden- 
ing effect of photometric errors, as derived from the test 
described above, then the implied intrinsic width of the 
red sequence is « 0.1 mag for all Zphnt 5, 1.2, c ons istent 
with the Zspec < 1-0 findings of Ruhland et al.l ()2008il . 

7. THE COLOR EVOLUTION 
OF MASSIVE RED GALAXIES 

Our next task is to quantify the color evolution of 
the massive red galaxy population. We have addressed 
this question based on the double Gaussian fits to the 
(CM*R-corrected) color distributions of > 10" Mg 
galaxies shown in Figure [6l In Figure [3 we plot the fit 
centers of the blue {lower points) and red {upper points) 
galaxy color distributions as a function of redshift; the 
errors on the locus of the massive red galaxy population 



shown in this Figure have been obtained by bootstrap 
resampling. 

From this plot, it is clear that the red galaxy pop- 
ulation as a whole has become progressively redder by 
~ 0.4 mag in [u — r) over the past 9 Gyr; the evo- 
lution in the blue cloud is similar. Making a linear 
fit to the (robust) Zphot < 1-1 measurements, we find 
A(u - r)corr = 2.57 - (0.44 ± 0.02) Zphot- (Note that 
this fit is constrained to match the z « point.) These 
results do not change significantly if we also include the 
point Zphot "^1.2, but it is clear that if we were to fit to 
the Zphot ^ 1-3 points, we would find slightly less strong 
evol ution. 

Kri ek et al.l (|2008D report a 3.3cr detection of a red se- 
quence in the spectrally derived {U — B) color distribu- 
tion of a mass-selected sample of 36 Zphot ^ 2 galaxies, 
12 of which lie in the ECDFS. The square at z = 2 in 
Figure [8] shows the approximate {u — r) color equivalent 
oj^their red sequence detection. While the Kriek et al] 
(j2008( ) point is slightly redder than an extrapolation of 
our linear fit, the two results agree rather well. 

In Figure [9l we compare the observed color evo- 
lution of the red galaxy population with nai've ex- 
pectations from passive evolution of synthetic spec- 
tra. For this purpose, we have us ed Pegase V2.0 
()Le Borgne fc Rocca-Volmerangd [2"002[ ) models with an 
initial metallicity oi Z = 0.004, and assuming a short 
burst of star formation (e-folding time of 100 Myr), be- 
ginning at Zform = 2, 3, 5, or 10 {bottom to top). 

The Zform <^ 3 tracks do an adequate job of describing 
the amount of observed evolution for Zphot < 1-3. How- 
ever, to get this level of agreement, the model (luminosity 
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weighted) metallicities at z = must be roughly solar 
{Z « 0.2), whereas the local mass-metallicity relation 
would suggest that these abu ndances should be sup er- 
solar by a factor of 3 or more ()Tremonti et all 120041 ). If 
we were to adopt Z — 0.2 initially, leading to final metal- 
licities that are super-solar by approximately 50 %, the 
model colors would become nearly 0.5 mag too red. The 
problem is even worse for SSP models: while solar metal- 
licities lead to colors which are too red by 0.25 mag, an 
acceptable description of the data is only possible assum- 
ing roughly half solar abundances. (See also Bell et al. 
2004.) 

Nevertheless, these simple models do provide an ac- 
ceptable description of the rate or amount of color evo- 
lution of the red galaxy population, A(m — r), if not the 
(u—r) colors per se. Using these extremely simple models 
to interpret the color evolution of the red galaxy popula- 
tion, the observations suggest that the bulk of the stars 
in massive red gal axies were fo r med at ^ ^ 3; includ- 
ing the point from iKriek etld] (j2008l ) would suggest a 
formation redshift > 5. 

This is not to say, however, that the observations are 
consistent with all massive red galaxies being formed at 
z > 5, or even 3: even among M* > 10"'^^ Mq galax- 
ies, there are simply not enough stars at z > 2 to 
build the z ~ red sequeii ce population (jFontana et al.l 
[20061: iConsehce et al.l[2007l see also Figure [T0|). Instead, 
what we see is that the colors of massive red galaxies 
are consistent with being dominated by ancient stars at 
all redshifts. This implies both extended star forma- 
tion histories among red galaxies (and/or their progen- 
itors), as well as a large spread in the times at which 
galaxies of a given (stellar) mass join the red popu- 
lation — a kind of long migration of gala xies, occur- 
ring over many Gvr (cf. iBrown et al.l[200il . A proper 
description of red sequence evolution would therefore 
have to account for, arnong o ther things, progenitor bias 
(|Van Dokkum fc Franx|[2001h : the continua l skewing of 
the population by new ad ditions (see also iFaber et ahl 
I2007t iRuhland et al]l2008f l. This is beyond the scope of 
this work. 



8. THE RISE OF RED GALAXIES OVER Zphot < 2 

In this section, we turn our attention to the third of 
our basic results. Whereas our focus until now has been 
on the properties of the red galaxy population as a whole, 
we now look at how the number of red galaxies within 
the total massive galaxy population — ie. the red galaxy 
fraction — evolves with time. 



8.1. Defining a Red Galaxy Selection Criterion 

Our inability to robustly distinguish separate red and 
blue galaxy populations on the basis of the observed color 
distributions for Zphot ^1-2 forces us to devise some al- 
ternate means of separating 'red' galaxies from the gen- 
eral field population. 

We have already seen that the color evolution of the 
red galaxy population is roughly consistent with ancient 
stars at all redshifts (Figure ^ . Our simple solution is 
therefore to use the predicted color evolution for a pas- 
sively evolving stellar population formed at high redshift 
to define a redshift-dependent 'red' selection criterion. 



VIZ.: 

{u-r)> 2.57 + 0.24 x logio(M,/10"MQ) 

+ 5(^phot)-0.25 , (2) 

where 5(zphot) is the (u — r) color evolution of a SSP with 
Zfoim = 5, as shown in Figure [HI This selection limit is 
shown as the dashed lines in each of Figures [5l [6l and [51 

How does this definition of 'red' relate to things like 
membership of the red sequence, star formation rate 
and/or history, etc.? As we remarked in the first para- 
graph, a galaxy's optical color is a reflection of its mean 
stellar age, modulo the complicating factors of metallicity 
and dust extinction. In addition to 'red and dead' galax- 
ies, therefore, simply selecting 'red' galaxies can poten- 
tially catch a significant number of star forming galaxies 
with high dust obscuration. In other words, while all 
passive galaxies are red, not all red galaxies are passive. 

In this sense, it is not unreasonable to interpret the 
redshift evolution of the number and fraction of red, mas- 
sive galaxies as placing an upper limit on the numbers 
of 'fully formed' (in the sense that they have essentially 
completed their star formation and/or assembly) massive 
galaxies. These results can thus be used to constrain the 
epoch at which the star formation quenching mechanism 
operates. 

8.2. The Number Density Evolution of Massive 
Galaxies for Zphot < 1-8 

Figure [TOj shows the evolving number density of > 
10^^ Mq galaxies for Zphot < 1-8. As before, the 
z « point comes from our analysis of the low-z sam- 
ple discussed in Appendix [XI the histograms are for the 
main MUSYC ECDFS sample. For the z > galax- 
ies, since we have used bins of equal comoving volume, 
the observed numbers (right axes) can be directly re- 
lated to a comoving number density (left axes), mod- 
ulo uncertainties in the cosmological model. The black 
point/histograms refer to the total > 10^^ Mq pop- 
ulation; the red point/histograms refer to red galaxies 
only; the dotted lines show where our results are signifi- 
cantly affected by incompleteness. 

The error bars on the z ^ histograms include 
the estimated measurement uncertainty due t o field -to- 
field variation, derived as in lSomerville et all ()2004l ). but 
modified for cuboid rather than spherical volumes (R 
Somcrville, priv. comm.). For any single measurement, 
this is the dominant source of uncertainty: typically ~ 30 
%, as compared to random photometric errors, which are 
at the ~ 10 % level. Note, however, that each of the 
•Zphot ^ 1 bins contains i ts own 'spike' in the Zgpec dis- 
tribution (Vanzellaet J][200l); the 0.6 < Zphot < 0.8 
bin contains two, and the 0.8 < Zphot < 1-0 bin none. 
Further, note that for Zphot ^ 1, our redshift errors are 
comparable to the size of the bins themselves; in this 
sense, the densities in neighboring bins are correlated, 
and the variations due to large scale structure are to a 
certain extent masked. 

The observed evolution in the number density of mas- 
sive galaxies is moderate: the number density of A/* > 
10" Mq galaxies at 1.5 < Zphot < 1-8 is 52 ± 8 % of 
the z « value. We see a very similar trend in the 
mass density, although a handful of galaxies with inferred 
Af* ^ 10^^ Mq make this measurement considerably 
noisier. 
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Fig. 10. — The rise of massive, red galaxies over cosmic time — Lower panel: evolution in the number density of all galaxies (black 
histograms) and of red galaxies (red histograms) with M, > 10^^ Mp in th e MUSYC ECDFS catalog; The error bars shown reflect 
the expected field— to-field variation, calculated as per ISomerville et alj II2004I) . Upper panel: evolution in the red galaxy fraction among 
M* > 10^^ M0 galaxies, as a function of photometric redshift. In this panel, the error bars have been derived by bootstrap resampling. 
In both panels, the dotted histograms show where we are significantly affected by incompleteness. The 2 « points (circles) are derived 
from the low-z sample, discussed in Appendix |A] The smooth curves show our fits to the observations; the shaded regions show the la 
uncertainties in the fits: we find 7tot = —0.52 it 0.12, 7rcd = —1.70 ± 0.14, and 7frac = —1.17 ± 0.18. 



In order to quantify the observed evolution, we have 
made a parametric fit to our measured number densities 
of the form: 

ntot(;2) = no (l + z^-' . (3) 

In practice, since the dominant source of uncertainty in 
any single point is from field-to-field variance, we per- 
form a linear fit in logn-log(l + z) space, weighting each 
point according to its uncertainty as shown in Figure 
[TUl Further, we constrain the fits to pass through the 
z « point, effectively eliminating uq as a free param- 
eter. In this way, we find 7tot — —0.52 ± 0.12. The fit 
itself is shown in Figure [10] as the smooth black curve; 



the shaded region shows the Icr uncertainty in the fit. 

8.3. Evolution in the Red Galaxy Fraction from 
•2^phot — 1-8 to the Present Day 

Wc now turn our attention to the red galaxy popula- 
tion. From Figure [TOl it is clear that the observed evo- 
lution is much stronger for red galaxies than it is for the 
total population: the number density of red galaxies at 
1-5 < Zphot < 1-8 is 18 ± 3 % of thc'z « value. Mak- 
ing a fit to the red galaxy number densities to quantify 
this evolution, in analogy to the previous section, we find 
7iod = — 1.60± 0.14 {smooth red curve). 
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A complementary way of characterizing the rise of mas- 
sive red galaxies is to look at the evolution of the red 
galaxy fraction. There are several advantages to focus- 
ing on the red galaxy fraction, rather than the number 
density of red galaxies per se. First and foremost, the 
uncertainty in the red fraction due to large scale struc- 
ture and field-to-field variation should be considerably 
smalle r than those in the number density, especially at 
z > 1 (|Cooper et al.ll2007[ ). 

We show the red fraction as a function of photometric 
redshift in the upper panel of Figure [TOl Fitting these 
results (smooth red curve) with the same functional form 
as in Equation [3l we find 7frac — —1-06 ± 0.16. Taken 
together, the results encapsulated in Figure [TO] suggest 
that < 20 % of > 10" M© galaxies in the local 
universe were already on the red sequence by Zphot ~ 1.6 
(9.5 Gyr ago). By the same token, approximately 50 % 
of these galaxies only (re-)joined the red sequence after 
z = 0.5 (5.0 Gyr ago). 

8.4. The Importance of the z w Comparison Point 

It is clear from Figure \W\ that almost all the signal 
for z ^ evolution in the red fraction comes fro m the 
z w comparison point (see also lBorch et al.|[2006f) . Fit- 
ting to the z ^ data alone, we find /red = (0.53 ± 
0.05) (1 -I- Zphot)"°-^^=^°-^^; that is, less than a 2a sig- 
nal of evolution. The same is a lso true for the nu mber 
density measurements (see also iBorch et al.l [2006) ; the 
reasons for this are not just the relatively mild evolu- 
tion, but also to the relative size of the error bars on 
the low- and high-redshift measurements. Fitting only 
to the z points, we find 7tot — — 1.55 ± 0.84 and 
7rcd ~ —1.77 ± 1.84; these fits 'overpredict' the z « 
number densities by factors of 2.2 and 1.3, respectively. 

In this sense, then, rather than quantifying the abso- 
lute evolution in the z ^ population, we are perform- 
ing a difference measurement between the situations at 
z « and z 3> 0. For this reason, it is imperative that 
care is taken when deriving the z w comparison values 
to ensure that the low- and high-redshift samples have 
been analyzed in a uniform way (Appendix [X)) . 

In com p arison to more soph i sticat ed analyses by 
iBeh et al.l (2003) and iGole et al.l (l200l . both of which 
combine the 2MASS and SDSS datasets, our z « num- 
ber densities are approximately 10 % higher, and 15 % 
lower, respectively. Adopting these values in place of our 
own deterrninati ons, we find 7tot — —0.72 ±0.12 using the 
IBeh et al.l (I2003D mas s function, and 7tot = -0.42 ± 0.12 
usmg the IGole et al.l (12001.') mass function, changes of 
—0.2 and -1-0.1 with respect to our default analysis. 
These rather large discrepancies underline the impor- 
tance of uniformity in the analysis of high- and low- 
redshift galaxies. 

8.5. Comparison with Other Works 

In Figure [Til we compare our results to a selec- 
tion of the steadily growing number of similar mea- 
surer nents; we show res ults from: the GOMBO-17 sur- 
vey fBorc h et al.l [2006h . the GOODS- MUSIG catalog 
jFontana et al the K20 surve y (jFontana et al.l 

[2004I). the VVDS (iPozzetti et al. ''2007). the DEEP2 sur- 
vey (iBundv et all 120061) . and^UNIGS (Drorv ,20041) . 
(Note that all of these results have come within the 



lookback time (Gyr) 
6 7" 



2 0.0008 




1.75 2.00 2.25 



photometric redshift, z. 



phot 

Fig. 11. — Comparison with other worlcs — The different symbols 
correspond to different authors and surveys as marlsed: the filled 
symbols refer to other authors' analyses of SDSS data; references 
for other z 3> surveys are given in the main text. As in Figure 
[TOl the smooth curves show our fits to the MUSYC ECDFS data; 
the shaded regions show the Icr errors on the fits. Excepting the 
MUSYC ECDFS data points, all points have been derived by inte- 
grating up Schechter function fits to the observed mass functions, 
in redshift bins. Apart from the MUNICS and the COMBO-17 
ECDFS data points, there is good agreement between the many 
different groups' results; the cause for the discrepancy between the 
COMBO-17 and MUSYC results in the ECDFS is investigated in 
Appendix |B] 



last five years.) In all cases, the points shown in Fig- 
ure [TT] have been obtained by integrating up a fit mass 
function, taking into account different choices of cosmol- 
ogy and IMF. We note that the strong redshift spike at 
Zphot ^ 0.7 is also present in the GOODS-MUSIG results, 
which are based on the GDFS. 

The comparison with the GOMBO-17 results deserves 
some further comment. While our results agree reason- 
ably well with the combined results from three of the 
GOMBO-17 fields {large diamonds), there is significant 
disagreement between our results and the GOMBO-17 
results in the EGDFS (small diamonds). In Appendix [Bl 
we present a detailed comparison between the MUSYG 
and GOMBO-17 datasets, and demonstrate that the 
discrepancy between the two results is due to system- 
atic calibration errors in the GOMBO-17 photometry, 
rather than differences in our analyses. We note in pass- 
ing that these calibration errors have only a small ef- 
fect on the GOMBO-17 photometric redshifts, and af- 
fect the ECDFS only; not the other GOMBO-17 fields 
(|Wolf et al.ll200l) . 

Particularly given the significant uncertainties in these 
measurements, and with the possible exception of the 
MUNIGS results, the degree of agreement between these 
surveys is impressive. However, in all of these surveys, 
field-to-field variance is a — if not the — major source of 
uncertainty. More and larger fields are necessary to lock 
down the growth rate of massive galaxies. 

9. QUANTIFYING POTENTIAL SYSTEMATIC 
ERRORS 
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Having now described our experiment in full, in this 
section we describe a wide array of sensitivity analyses, 
in order to determine how sensitively our results depend 
on specific choices in our experimental design. The basic 
idea is to vary individual aspects of our analysis, and see 
what effects these changes have on the results presented 
in iJSl in particular, we will focus on how the best-fit val- 
ues of 7tot and 7frac depend on our experimental design. 
The results of many of these tests are summarized in Ta- 
bleslHandO Note that we have used the same red galaxy 
selection criterion (Equation [2]) throughout all of these 
tests. 

We will consider, in turn, how sensitive or robust our 
results are to basic photometric calibration errors fi i9.ip . 
the methods used for the photometry itself f ii9.2p . in- 
completeness ( TO.3p . and variations on our photometric 
redshift computation ('i i9.4p . We then discuss possible 
systematic effects arising from our method for estimat- 
ing stellar masses in TO 51 and other concerns in TO. 61 
Our findings in the section are summarized in TO. 71 

9.1. Photometric Calibration 

How sensitive or robust are our results to errors in the 
basic photometric calibration? And how accurate are the 
calibrations of each of the different bands, relative to one 
another, and in an absolute sense? We address these two 
questions in this section. 

9.1.1. The Effect of Perturbing Individual Bands 

In order to gauge the effect of calibration errors on our 
final results, we have perturbed the photometry in each 
band in turn by ±0.05 mag, and repeated our full anal- 
ysis, from the derivation of photometric redshifts and 
stellar masses to fitting the 7s. Roughly speaking, these 
perturbations can affect our results in two separate ways: 
either through direct changes in the SEDs themselves, 
or indirectly, through changes in the derived photomet- 
ric redshifts, and so the transformation from observed to 
restframe quantities. In order to disentangle these two ef- 
fects, we have also repeated our analysis with these ±0.05 
mag shifts, but while holding the photometric redshifts 
fixed to their default values. The results of these tests 
are summarized in Table [2l 

For the U BV bands, the effect of changing the photo- 
metric calibration on our main results is driven almost 
exclusively by changes in the photometric redshifts. This 
is simply because our stellar mass estimates typically do 
not depend directly on the observed U BV photometry: 
the restframe B has already shifted past the observed 
V hy z ~ 0.25. (Recall that we use the B — V color 
to infer M^/Ly, and thus M*.) Similarly, the effect that 
changing the reddest H and K bands has on 7tot is small: 
± < 0.02 for both. The situation changes for the RIz'J 
bands, where the direct effect comes to dominate over 
the indirect effect from changes in the redshifts. 

In terms of the SED photometry, our results are most 
sensitive to errors in the R and J band photometric cal- 
ibrations: a ±0.05 mag shift in the R oi J zeropoint 
changes the value of 7tot by =f0-08 or ±0.14, respectively. 
For the R band, the effect on restframe properties is fo- 
cused on the Zphot ~ 0.7 bin, where we happen to mea- 
sure the highest density; this is also the z ^ bin that 
has the single greatest effect on our measurement of 7tot- 
The critical importance of the J band stems from the 



fact that it plays a role in the derivation of restframe V 
photometry, and thus M*, for all Zphot ^ 1; ie-, roughly 
two thirds of our surveyed volume. It seems plausible 
that our sensitivity in this regard might be substantially 
reduced if we were to use a different method for estimat- 
ing stellar masses (see also TO.5p . 

While the SED modeling results are not particularly 
sensitive to the K band calibration, it still plays a critical 
role in normalizing each SED through the total K flux 
measurement. Although changing the K calibration by 
±0.05 mag effectively rescales the stellar masses by just 
=F0.02 dex, this can change the inferred number densities 
by as much as =FlO %■ In terms of 7tot, the effect is 
TO.07. 

9.1.2. Testing the MUSYC ECDFS Photometric 
Calibration 

We now turn our attention to identifying and quanti- 
fying potential calibration errors in the MUSYC ECDFS 
dataset. We will then be able to use this information to 
estimate the extent to which our results may actually be 
affected by such errors. 

To address this question, we have used EAZY to fit 
main sequence stellar sp ectra from the BPGS stellar atlas 
(|Gunn fc Strvkei|[l983( ). fixing Zphot = 0, to the observed 
photometry for stars. Assuming that whatever calibra- 
tion errors do exist do not affect the choice for the best 
fit template spectrum, we can then interpret the mean 
difference between the best fit and the observed photom- 
etry as being the product of calibration errors. We note 
that since this test is concerned exclusively with SED 
shapes, it cannot comment on the absolute calibration of 
any given band; instead it assesses the relative- or cross- 
calibration of the ten band photometry that makes up 
the MUSYC ECDFS dataset. 

The offsets we derive in this way are given in Col. (10) 
of Table [2] The most notable offsets derived in this way 
are -0.048 in (CZ-C/ss), 0.033 in {I-z'), 0.064 in {J-H), 
and 0.032 in {J — K). These offsets give an indication of 
the plausible size of any calibration errors or inconsisten- 
cies in the MUSYC SED photometry. If we recalibrate 
our photometry to eliminate these apparent offsets (hold- 
ing the K band fixed), and repeat our full analysis, we 
find that the best fit values of 7tot , 7rcdj and 7frac change 
by +0.00, +0.02, and +0.04, respectively. 

We have also checked the absolute calibration of 
each band in comp arison to the FIREWORKS catalog 
(jWuvts et al.ll2008l ) in the region of overlap (Paper I) . Us- 
ing (B — z')-{z' — K) selected stars, we have constructed 
empirical 'color transforms' between FIREWORKS and 
MUSYC filters, as a function of HST color. Comparing 
these diagrams to pr edictions derived from the BPGS 
stellar spectral atlas ()Gunn fc StrvkejfTgsl . we treat 
any offset between the predicted and observed stellar se- 
quences as a calibration error in the MUSYC photome- 
try. (This is the same method we used to identify the 
disagreement between the COMBO-17 and MUSYC cal- 
ibrations discussed in Appendix IbI) 

The offsets we have derived in this way are given in Col. 
(11) of Table [2] The biggest offsets derived in this way 
are « 0.05 mag in (B-F435W) and (/ - F850LP). The 
offset in (i^-i^isAAc) is —0.017 mag; for the crucial J 
band, the offset is just +0.015 mag. In agreement with 
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TABLE 2 

Photometric Calibration Sensitivity Tests 





^phot 


fixed 






Recalculating Zpi^ot 






Photometric offset with respect to: 




change 


in 7tot 


change 


in 7tot 


change 


in 7red 


change 


in 7frac 


stellar SEDs 


FIREWORKS 




-fO.05 mag 


—0.05 mag 


+0.05 mag 


—0.05 mag 


+0.05 mag 


—0.05 mag 


+0.05 mag 


—0.05 mag 








(2) (3) 


(4) (5) 


(6) (7) 


(8) (9) 


(10) 


(11) 


TJ 
U 


+0.00 


+0.00 


+0.01 


-0.01 


+0.02 


+0.02 


+0.001 


+0.041 


—0.004 


+0.013 


U3a 


+0.00 


+0.00 


+0.00 


+0.01 


+0.01 


+0.02 


-0.002 


+0.009 


-0.051 


-0.020 


B 


+0.00 


+0.00 


+0.00 


+0.00 


+0.03 


+0.01 


+0.061 


-0.006 


-0.017 


+0.050 


V 


+0.01 


+0.00 


-0.02 


+0.00 


+0.04 


-0.10 


+0.037 


-0.068 


-0.006 


+0.038 


R 


+0.05 


-0.03 


+0.08 


-0.08 


+0.06 


-0.05 


+0.071 


-0.053 


+0.017 


+0.016 


I 


+0.04 


-0.03 


+0.01 


-0.02 


+0.03 


-0.03 


-0.051 


+0.066 


+0.023 


+0.055 


z' 


-0.02 


+0.01 


-0.03 


+0.02 


-0.01 


-0.01 


-0.011 


-0.024 


-0.011 


-0.004 


J 


-0.10 


+0.09 


-0.13 


+0.10 


-0.14 


+0.11 


+0.022 


-0.033 


+0.032 


+0.015 


H 


-0.03 


+0.03 


-0.02 


+0.01 


-0.02 


+0.05 


-0.009 


+0.039 


-0.032 


-0.012 




+0.00 


+0.02 


-0.02 


+0.01 


-0.07 


+0.08 


-0.078 


+0.096 






Ktot 


-0.07 


+0.07 


-0.07 


+0.07 


-0.03 


+0.03 


+0.057 


-0.025 




-0.017 



Note. — For each band (Col. 1), Col.s (2) — (9) give how our parametric fits to the 2phot ^ 2 evolution of the Af» > 10^^ Mq popuation 
changes with ±0.05 mag perturbations to individual bands' zeropoints: Col.s (2) and (3) show the effect on the total number density 
evolution of massive galaxies due to the change in the photometry only. Col.s (4) and (5) show the same information, but including the 
effect of changes in the photometric redshifts that come with changing the calibration; similarly, Col.s (6) and (7) show the effect on the 
number density measurements for red, massive galaxies, and Col.s (8) and (9) show the effect n the red galaxy frac tion. Col.s (1 0) and 
(11) give zeropoint offsets for each band based on stellar colors and comparison with the FIREWORKS catalogue l|Wuvts et al.l 120081 . 
respectively. 



the stellar SED fitting test, this analysis also finds an 
inconsistency between the I and z bands at the level of 
0.05 mag. If we recalibrate our photometry to match the 
FIREWORKS catalog, we find that our values of 7tot, 
7red, and 7frac changc by —0.03, —0.13, and —0.08 with 
respect to our default results. The sizes of these changes 
are in excellent agreement with the results of the previous 
section; more than half of these changes can be ascribed 
just to the 0.02 mag rescaling of total K magnitudes. 

We estimate that the systematic uncertainty in our 
main results due to photometric calibration errors is at 
the level of A7tot < 0.05 and A7ft.ac 0.10. 

9.2. Photometric Methods 

While we rely on SExtractor for our basic photometry, 
we have introduced three sophistications in our analysis. 
In this section, we investigate the effects that these three 
changes have on our results. 

9.2.1. Background Subtraction 

We have applied a correction to account for the ten- 
dency of SExtractor to overestimate the background 
level; for individual objects, this correction is typically 
on the order of —0.03 mag. To test the sensitivity of 
our results to background subtraction errors, we have 
repeated our analysis without applying this correction. 
Repeating our analysis relying on our own background 
estimation we find changes in 7tot and 7i.od of —0.06 and 
—0.05, respectively; the change in 7frac is just —0.02. 

9.2.2. Total Magnitudes 

Our total flux measurements are based on SExtrac- 
tor 's AUTO flux measurement. It is well known that 
the AUTO aperture misses a significa nt amount of flux, 
espec i ally for the faintest sources (iBertin fc Arnouta 
19961: iLabbe et al.[ l2003l: iBrown et al.ll2007D . Following 
Labbe et al.l (|2003f) . we partially redress this by applying 



a minimal correction to correct for light laying beyond 
the AUTO aperture, treating each object as if it were 
a point source. For individual sources, this correction 
amounts to as much as 20 %. Experiments with syn- 
thetic r^/^ sources placed in our own image suggest that 
for the specific example of a i^T = 22 elliptical galaxies 
with i?e = 0'.'4 (w 3 kpc at z = 1), this correction reduces 
the missed flux from 0.33 mag to 0.15 mag. 

In terms of the measured number densities, the effect of 
adopting this correction is 1 — 5 % for Zphot ^1-5, rising 
to 5 — 10 % for 1.5 < Zphot < 1-8. The use of uncorrected 
AUTO fluxes as total flux measurements thus produces a 
rather mild spurious evolutionary signal. Repeating our 
analysis without this correction, we find slightly steeper 
evolution in the number densities: the values of 7tot, 
7rcd change by —0.07 and —0.06, respectively. By com- 
parison, the red fraction measurements remain almost 
unchanged: 7frac changes by —0.01. 

9.2.3. SED Apertures and Colour Gradients 

When constructing multi-color SEDs, we have used the 
larger of SExtractor's ISO aperture (based on the I'.'O 
FWHM K mosaic) and a fixed 2'.'5 diameter aperture. 
This aperture flexibility is intended to guard against po- 
tential biases due to color gradients in individual galax- 
ies. This problem is presumably the most severe for 
the largest, relatively low redshift galaxies with signifi- 
cant bulge components, leading to overestimates in both 
{B - V) and M,. 

Consistent with this expectation, when repeating our 
analysis relying exclusively on the fixed 2'.' 5 aperture 
fluxes to construct SEDs, we find the measured num- 
ber densities of massive galaxies increases by 5 % for 
•Zphot < l-l- The increases in these low redshift bins 
brings them closer to the z ^ point, leading to a slight 
decrease in the measured evolution: 7tot and 7].od change 
by +0.05 and -1-0.13, respectively; 7trac changes by just 
+0.02. 
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We estimate that the systematic uncertainty in our 
results associated with our photometric methods are on 
the order of 7tot ^ 0.07 and 7fi.ac ~ 0.00. 

9.3. Correcting for Incompleteness 

In i|4l we have argued that we are approximately com- 
plete for M* > 10^^ M© galaxies for Zphot < 1-8. In this 
section, we examine the effects of incompleteness due to 
both the K < 22 detection threshold and the K S:N 
> 5 'analysis' selection, by deriving simple completeness 
corrections. 

9.3.1. Incompleteness Due to Low S:N 

In addition to incompleteness due to our K < 22 de- 
tection limit, discussed below, there is the concern of in- 
completeness due to the if S:N > 5 selection limit, which 
we have imposed to ensure against extremely poorly con- 
strained photometric redshifts. This cut affects 4.5 % of 
all K < 22 sources in the ECDFS catalog, with 43 % 
of those galaxies laying in the slightly shallower Eastern 
pointing. At a fixed magnitude, there is not a strong 
dependence of the fraction of S:N < 5 detections as a 
function oi J — K color or — with the caveat that this cut 
is designed to remove poorly constrained redshifts — as a 
function of Zphot- 

In order to assess the impact of this cut on our results, 
in this section we attempt to correct for this additional 
source of incompleteness. Our procedure is as follows: we 
determine the completeness fraction in the face of this se- 
lection, as function of total K magnitude, f{K); this can 
then be used to weight each retained galaxy according to 
its K magnitude, w{K) — l/f{K). Of galaxies in our 
main sample, 8 % of all galaxies, and 11 % of > 10^^ 
Mq galaxies at 1 < Zphot < 2 are given w{K) > 1.25. 

Recalculating our number density measurements using 
these weights, we find that the number density of mas- 
sive galaxies at 1.2 < Zphot < 1-8 changes by less than 5 
%; the same is true for red galaxies alone. Repeating our 
analysis with these corrections, we find 7tot = —0.49, 
7red = —1.56, and 7frac = —1.05, amounting to differ- 
ences of -1-0.03, -1-0.04, and -fO.Ol with respect to our 
fiducial analysis. We get similar results adopting a more 
stringent S:N > 10 criterion. To put this change into 
perspective, it is comparable to that due to uncertainty 
in the background subtraction for our basic photometry. 

9.3.2. Correcting for Undetected Sources 

We have corrected for incompleteness due to our K < 
22 detection limit as follows: taking the Af* > 10^^ Mq 
galaxy population in a given redshift bin, we then predict 
what the observed K magnitude for each galaxy would be 
if the galaxy were shifted through the next most distant 
redshift bin. In other words, we effectively 'K correct' 
each galaxy's observed K flux from its fiducial Zphot to 
higher redshifts, using the same machinery as for the 
interpolation of restframe fluxes. This makes it possible 
to determine the fraction (in terms of volume) of the 
next redshift bin over which a given galaxy would remain 
detectable; the overall completeness is then simply the 
average of this quantity for all galaxies in the bin. We 
have performed this correction for red and blue galaxies 
separately. 

Based on this analysis, we are indeed 100 % complete 
for Zphot < 1.25; for 1.7 < Zphot < 1-8, we are more 



than 80 % complete overall, and at least 75 % complete 
for red galaxies. This agrees reasonably well with our 
completeness estimates in gl For higher redshifts, our 
estimated completeness drops to 75 % for 1.9 < 2:phot < 
2.0, and to 65 % for 2.1 < Zphot < 2.2. The estimated 
incompleteness correction factor to the measured number 
densities is thus less than 1.5 for all 2;phot < 2.2. Using 
these completeness corrections to extend our analysis to 
2phot < 2.2, we find 7tot = -0.47 and 7frac = -0.90; 
changes of —0.05 and —0.16 with respect to our default 
values for Zphot < 1.8. 

Fitting to these incompleteness-corrected measure- 
ments for the number density of massive galaxies with 
•Zphot < 1-8 {ie., repeating our default analysis, but with 
a completeness correction), we find 7tot = —0.37, —1.42, 
and 7frac = —1.06. In comparison to our uncorrected 
results for Zphot < 1-8, these represent changes of —0.15, 
—0.18, and —0.00(3), respectively. However, we also note 
that if we were to assume that we are 100 % complete 
for Zphot < 1-4, these changes become —0.10, —0.12, and 
-f 0.01; that is, the uncertainties on these incompleteness 
corrections are comparable to in size to the corrections 
themselves. 

9.4. Photometric Redshifts 

The next major aspect of our analysis that we will ex- 
plore in detail is systematic effects associated with the 
derivation of photometric redshifts; we split this discus- 
sion into two parts. In the first part part f ^9.4.ip we 
explore how our results depend on the choice of tem- 
plates used in the Zphot calculation. Then, in the second 
(i j9.4.2[) . we investigate how our results depend on the ex- 
act method used for deriving photometric redshifts, by 
varying individual aspects of the EAZY algorithm. Some 
illustrative results from a selection of these sensitivity 
tests are given in Figure \T% 

9.4.1. Trialing Different Template Sets 
Using the \Fontana et al\ 1(200 A ) template set — 



Using a library of ^ 3000 sy i itheti c Pegase V2.0 
(Lc B orgne fc Rocca-Volmeran"gp ^ I2002D spectra de - 
scribed bv iFontana et all (|2006f) . iGrazian et all (|2006f ) 
obtain « 0.045 for K selected galaxies from their 
catalog of the GOODS-CDFS data. In fact, this library 
provides the parent catalog for EAZY's default template 
set (see Brammcr et al. 2008). If we use this template 
library in place of the EAZY default, and do not allow 
template combination, we find CTz — 0.039; the fraction 
of objects having Az/{1 + z) > 0.1, /o.i, is 16 % 
(c/. 12.4 % for our default redshifts). Adopting these 
redshifts in place of our default determinations, we find 
7tot = —0.55 and 7ti.ac — —1.00; differences of —0.03 and 
—0.06, respe ctively (see Figure [T^ . 

Using the \Rudnick et al\ \200A ) template set — In the 
past, our group ha s tended to d e termi ne photomet- 
ric redshifts as per iRudnick et al.l ()2003D . which con- 
siders linear combin ations among: four empirical tem - 
plate spectra from iColeman. Wu fc Weedmaiil (|1980D : 
two starburst spectra from iKinnev et al.l ()1996( ): and 
two lBruzual fc Chariot) (|2003l ) synthetic spectra for SSPs 
with ages of 10 Myr and 1 Gyr. If we use these templates 
with EAZY, we find ct^ = 0.055, with /o.i = 27 %. Look- 
ing at Figure [T^ these photometric redshift determina- 
tions suffer from a number of systematics whereby there 
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Fig. 12. — The effect of different photometri c redshift analyse s on our results — In rows, we show how our results would change {top to 
bottom): adopting the template set used by Rud nick et al.l (120031): adopting the template set described by ?; emulating hyperz in its default 
configuration; adopting the maximum likelihood {ie.y^ minimisation) photometric redshift value; neglecting the U band photometry. In 
each row, we show {left to right) the analogues of Figure [T] the third panel of Figure |5] and Figure [TOl 
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are clear preferred redshift solutions at, e.g., 2;phot ^ 0.4, 
as well as a larger systematic underestimate of the true 
redshift: in addition to the larger random error, we also 
find a systematic offset in Az/(1 + z) of —0.054. The 
effect this has on the Zphot ^ 0.8 number densities is 
strong, and is dominated by the shifting of many of the 
galaxies from the 0.6 < Zphot ^0.8 bin down into the 
0.2 < zphot < 0.6 bin. The net effect of this change 
on our final results is considerably less, however: using 
these redshifts, we find 7tot = —0.58 and 7frac = —1.11; 
changes of —0.06 and —0.05 with respect to our fiducial 
results. 

Emulating hyperz with empirical templates — As a fi- 
nal test for this section, we have also repeated our 
photometric redshift analysis emulating the popular hy- 
perz (Bolzpjiella, Miralles fc Pell o 2000.) code in its de- 
fault configuration; that is, minimization between the 
obse rved photometry and synthetic p hotometry using 
four I Coleman. Wu fc Weedmaiil ()1980D empirical tem- 
plate spectra, and making no allowance for dust extinc- 
tion beyond what is included in the empirical spectra. 
Using the spectroscopic redshift sample shown in Figure 
m we find ct^ = 0.060, and again a substantial systematic 
offset: Az/(1 + z) ^ -0.051. 

One very striking consequence of using so few template 
spectra is a very pronounced biasing of galaxy restframe 
colors, corresponding to the colors of the template spec- 
tra themselves (see Figure [T^ . This of course has a very 
strong effect on our stellar mass estimates. For this rea- 
son, it is critical when relying on photometric redshifts 
to derive restframe properties to use a template basis set 
that spans the full range of galaxy colors. 

9.4.2. Variations in our Default Photometric Redshift 

Calculation 

Adopting the best fit redshift — By default, EAZY as- 
signs each object a redshift by taking a probability 
weighted integral over the full redshift grid. However, 
it also outputs the single most likely redshift, deter- 
mined by minimization. For these redshifts, we find 



cr^ = 0.035, and /o.i = 14.7 %. For 1 < 



-spec 



< 2, 



however, we find CTz — 0.086. We also see evidence for 
rather strong systematic effects in the photometric red- 
shifts, such that there are preferred photometric redshift 
solutions, for example at Zphot ~ 0.8 and 1.2 (see Fig- 
ure I12[) , corresponding to the points where the optical 
breaks fall between the observed filters. With these red- 
shifts, the values of 7tot, 7rod, and 7frac change by —0.11, 
—0.05, and —0.14, respectively. 

No Luminosity Prior — Our default analysis makes use 
of a Bayesian K band luminosity prior. If we do not 
include this feature, we find az = 0.038 and fo.i — 13.5 
%. The best fit values of 7tot, 7rcd, and 7frac change by 
—0.11, —0.38, and —0.54, respectively. These changes 
are large; however, we note that these results are not 
internally consistent, in that it is not true that 7frac ~ 
7rcd/7tot- Further, these results are not consistent with 
the test we describe in iJlOl 

No Template Error Function — A novel feature of the 
EAZY code is the inclusion of a template error func- 
tion; ie. a systematic error, as a function of restrame 
wavelength, designed to down-weight those parts of the 
spectra like the restframe near-UV where galaxies show 
an greater intrinsic variation. If we do not include this 



feature, we find cr^ = 0.034, and /o.i = 15.3 %; the best 
fit values of 7tot, 7frac, and 7frac change by -0.04, -0.14, 
and —0.13, respectively. 

Two Template Combination — By default, we allow 
non-negative combinations of all six EAZY templates 
when fitting to the observed SEDs to determine the red- 
shift. If we instead allow combinations of only two (but 
any two) of the template spectra, we find ct^ = 0.038 and 
/o.i — 14.4 %. The fit values of 7tot and 7frac change by 
+0.02 and -0.02, respectively. 

'Best Effort' Photometric Redshifts — It is well known 
that the WFI U filter suffers from a red leak beyond 8000 
A. We find that we get our best Zphot-^spec agreement if 
we ignore the U band photometry in the computation 
of photometric redshifts. For the spectroscopic sample 
shown in Figures [1] and [21 we find — 0.035 and /o.i = 
12.0 %; this translates into an uncertainty in the stellar 
mass estimates due to redshift uncertainties of just 0.1 
dex. This only modest improvement in our photometric 
redshift estimates changes the fit values of 7tot and 7trac 
change by —0.01 and -1-0.05, respectively. 

These tests show that our results do depend rather 
sensitively on the details of the photometric redshift cal- 
culation, typically at the level of A7 < 0.15. In par- 
ticular, the inclusion of a luminosity prior is crucial in 
shaping our results. However, we note that in most 
cases there are objective reasons to favor our default red- 
shifts and/or results. Similarly, the use of different tem- 
plate sets changes our results at the level of A7 w 0.06. 
The effect of minor refinements to our photometric red- 
shift calculation (two template combination versus un- 
restricted; the inclusion/omission of the U band data) 
is small: A7 < 0.05. Note that, in any case, changing 
the photometric redshift calculation cannot change our 
qualitative results: a gradual increase in the red galaxy 
fraction with time. 

9.5. Stellar Mass Estimates 

We have relied on a rather simple method for esti- 
mating stellar mass-to-light ratios, based on restframe 
(B — V) colors. Certainly this method does not make use 
of the full amount of information available in the full, 
ten-band SED; on the other hand, it does ensure that 
the same information is used for all galaxies, irrespective 
of redshift. 

Nevertheless, we have trialed reanalyzing our data us- 
ing standard population synthesis SED fitting techniques 
to estimate stellar masses for our main sample, in order 
to assess the potential scale of biases arising from this as- 
pect of our experimental design. Specifically, we have fit 
iBruzual fc Chariot! (|2003l ) synthetic spectra to the ob- 
served photometry, with redshifts fixed to the fiducial 
•Zphot, using the hyperzspec utility, which is a part of the 
hyperz vl.2 release package (M Bolzonella, priv. comm.). 

While the random scatter in comparison to the color- 
derived stellar mass estimates is 0.18 dex, after correct- 
ing for IMF, cosmology, etc., the SED-fit masses are 
systematically lower by 0.3 dex. This is also true for 
both the FIREWORKS (iWuvts et al.|[200l or GOODS- 
MUSIC (IGrazian et al.l I2OO6D catalogs of the CDFS- 
GOODS data. Besides this offset, we do not see evi- 
dence for strong evolution in the normalization of the 
color relation in comparison to the SED fit masses (but 
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TABLE 3 

Summary of Sensitivity Test Results 





1.5 < Zphot < 1-8 




2phot < 1-8 




z^2 


Test description 


n.tot/"0 "rod/^O 


7tot 


Trcd 


"Tfrac 


"tot/'^O "■rcd/'^O 


(1) 


(2) (3) 


(4) 


(5) 


(6) 


(7) (8) 



Fiducial Results 



Default analysis (zr,\^[.f < 1.8) 


0.52 


± 0.08 


0.18 


± 0.03 


-0 


52 


± 





12 


-1 


60 


± 





14 


-1 


06 


± 





16 





57 





15 


Correcting for incompleteness (^phot < 2.2) 





57 





20 


-0 


47 


± 





09 


-1 


39 


± 





10 


-0 


90 


± 





11 





61 





19 


Photometric Calibration (^9.1.21 
















































Recalibration to match stellar SEDs 





54 





18 


-0 


52 


± 





12 


-1 


58 


± 





14 


-1 


02 


± 





17 





58 





16 


Recalibration to match FIREWORKS'*) 





49 





15 


-0 


55 


± 





12 


-1 


73 


± 





14 


-1 


14 


± 





15 





56 





13 


Photometric Methods (i\9.2l 
















































No correction for background oversubtraction 





49 





16 


-0 


58 


± 





12 


-1 


65 


± 





14 


-1 


08 


± 





17 





54 





15 


No correction for missed flux 





49 





16 


-0 


59 


± 





12 


-1 


66 


± 





14 


-1 


07 


± 





17 





53 





14 


Fixed 2'.' 5 apertures for SEDs 





52 





17 


-0 


47 


± 





12 


-1 


47 


± 





14 


-1 


04 


± 





14 





61 





18 




Incompleteness Corrections ( ^9.311 
















































Correcting for incompleteness (zphot < 1-8) 





57 





20 


-0 


37 


± 





12 


-1 


42 


± 





13 


-1 


06 


± 





15 





67 





19 


Correcting for K S:N < 5 incompleteness 





54 





18 


-0 


49 


± 





12 


-1 


56 







14 


-1 


05 


± 





16 





59 





16 


Usina Other z^h^t Template Sets (^9A.ll 
















































Using Rudnick et al. f2003) templates 





52 





14 


-0 


58 


± 





12 


-1 


57 


± 





14 


-1 


11 


± 





15 





54 





16 


Using Fontana et al. (2006) templates'*' 





53 





21 


-0 


55 


± 





12 


-1 


57 


± 





14 


-1 


00 


± 





09 





56 





16 


Emulating hyperz with CWW (1980) templates 





48 





24 


-0 


39 


± 





14 


-0 


83 


± 





15 


-0 


69 


± 





15 





66 





35 


Variations on our Default z^hr^t Analysis (^9A.2t 
















































Adopting the best fit redshift 





55 





22 


-0 


63 


± 





13 


-1 


65 


± 





14 


-0 


92 


± 





14 





51 





15 


No K luminosity prior 





41 





09 


-0 


63 


± 





12 


-1 


98 


± 





15 


-1 


59 


± 





22 





51 





10 


No template error function 





53 





15 


-0 


56 


± 





12 


-1 


74 


± 





14 


-1 


19 


± 





20 





55 





13 


Two template combinations 





54 





16 


-0 


50 


± 





12 


-1 


56 


± 





14 


-1 


08 


± 





15 





59 





16 


Neglecting U data ('best effort' ^phots) 





46 





19 


-0 


53 


± 





12 


-1 


59 


± 





14 


-1 


01 


± 





16 





57 





16 



Other Variations ([[93] and ^93t 



SED-fit stellar masses 





48 





19 


-0 


49 


± 





12 


-1 


62 


± 





14 


-1 


11 


± 





14 





59 





15 


Including spectroscopic redshifts 





51 





18 


-0 


54 


± 





12 


-1 


60 


± 





14 


-1 


06 


± 





14 





56 





15 


Excluding X-ray detections 





50 





18 


-0 


59 


± 





12 


-1 


63 


± 





14 


-1 


03 


± 





15 





53 





15 


GOODS area only 





47 





22 


-0 


63 


± 





17 


-1 


37 


± 





24 


-1 


00 


± 





25 





51 





20 



Note. — Each row describes a different sensitivity test, designed to determine how and how much our main results depend on our 
choice of experimental design. These tests are described in {[9] Each row gives an identifier (Col. 1), a different test; references are given 
for a description of each. For the 'Dissecting our Default Zphot Analysis tests, additional Zphot information is given in ij9.4.2l For each 
test, we give the measured number density of A/» > 10^^ Mq galaxies with 1.5 < Zphot < 1-8: these values are given for the total massive 
galaxy samples (Col. 2), and the red subpopulation (Col. 3), and are normalised relative to the z value for all galaxies. We also give 
the best fit values for the total 0.2 < Zphot < 1-8 evolution in the number density of all massive galaxies (Col. 4), in the number density of 
red, massive galaxies (Col. 5), and in the red fraction among massive galaxies (Col. 6). Finally, we give values for the number density of 
massive galaxies (Col. 7) and of red, massive galaxies (Col. 8), obtained by extrapolating the fits to z = 2. Dominant effects are marked 
with an asterisk; a summary of ourinterpretations of these results can be found in i|9.7l 



see lLin et al.ll2007l ). Since we have not been able to iden- 
tify the source of this offset, we have simply corrected for 
it. Note, however, that we have not refit the z « SEDs 
for this test. 

Using these SED-fit stellar masses, we find 7tot = 
—0.49 7rcd = —1-62, and 7trac = — l-H- These repre- 
sent differences of just -1-0.03, —0.02, and —0.05 with 
respect to our fiducial values. Similarly, usin g the pre- 
script ion for M/L as a functio n of (B—V) from Bell et ahl 
(|2003( ) rather than that from lBell fc de Jond (|2001h . we 
find 7tot = -0.41, 7i.cd = -1.55, and 7fi.ac = -1.12; 
changes of -1-0.11, -1-0.05, and —0.06, respectively. 

This suggests that our results are not grossly effected 
by our choice of method for estimating stellar masses, 



and especially not by the use of color- as opposed to 
SED-derived stellar masses. 

9.6. Other Potential Systematic Effects 

Large Scale Structure and Field-to- Field Variance — 
One interesting test of the effect of field-to-field vari- 
ance is to restrict our attention only to the MUSYC 
ECDFS coverage of the GOODS-CDFS region. For this 
143 □' sub-field, we find 7tot = —0.62, 7i.cd = —1.37, and 
7f,^^ = -1.00, differences of -0.10, -H0.23, and +0.06 
with respect to the full 816 □' ECDFS. In comparison 
to the full ECDFS, the dearth of galaxies at Zphot ~ 0.8 
is even more pronounced, and seems to extend over the 
range 0.8 < Zphot < 1-2; this appears to be the main 



The Rise of Massive Red Galaxies 



23 



driver of the strong change in 7i.od ■ Conversely, if we ex- 
clude the significantly underdense GOODS area, we find 
7tot = -0.49, 7rcd = -1.60, and 7frac = -1.10; differ- 
ences of -f0.03, —0.00, and —0.04, respectively. Again, 
this analysis underscores the fact that measurements of 
the red fraction are more robust against the effects of 
large scale structure than those of the number densities 
themselves. 

The z « value — In i i8.41 we have shown that the 
z w comparison point is critical in providing most of 
the signal for z ^ evolution in the massive galaxy pop- 
ulation. We have also considered our results vary if we 
change the way we treat the z w point. In our default 
analysis, we constrain the fits so that they pass through 
the z « points, effectively using the z « point to nor- 
malize the high-redshift measurements. Since the (Pois- 
son statistical) errors on the z ~ points are just a few 
percent, our results do not change significantly if we in- 
clude the z w points in the fits. 

For our default analysis, we have approximately ac- 
counted for the effect of photometric redshift errors, as 
they apply to the z ^ sample, on the z ~ mea- 
surements; we have done this by randomly perturbing 
the masses and colors of z « galaxies using the typ- 
ical uncertainties for z ^ galaxies, due to the use of 
photometric redshifts, and given in Figure [21 Systematic 
and random errors in the number densities both of all 
and of red galaxies are on the order of 5 %. If we do 
not account for the Eddington bias in this way, we find 
that 7tot — —0.46 and 7frac — —1.16, changes of 0.06 and 
-0.10, respectively. 

Correcting for Galactic Dust Extinction — Note that we 
have not included specific corrections for Galactic fore- 
ground dust extinction. The CDFS was specifically cho- 
sen for its very low Galactic gas and dust column den- 
sity; the suggested corrections for the optical bands are 
< 0.05 mag. These corrections are typically as large or 
larger than the uncertainties on the photometric zero- 
points themselves; it is for this reason that we have cho- 
sen not to apply these corrections. If we were to include 
these corrections, however, our derived values of 7tot and 
7frac change by —0.03 and -f 0.04, respectively. 

Including Spectroscopic Redshift Determinations — 
From among the thousands of spectroscopic redshift de- 
terminations that are publicly available in the ECDFS, 
we have robust ZgpccS for just over 20 % (1669/7840) of 
galaxies in our main sample, and for 20 % (269/1297) 
of those with 0.2 < Zphot < 2.0 and > lO" Mq. 
Repeating our analysis with this additional information 
included, our results do not change significantly: we find 
7tot = —0.50 and 7rcd = —1.70, differences of just -1-0.02 
and -1-0.00, respectively. Similarly, we show in Appendix 
[B]that our z < 0.8 results do not change by more than 
a few percent if we adopt photometric redshifts from 
COMBO-17 (Wolf et al. 2004; = 0.035). 

Excluding X-ray-Selected Galaxies — Note that we have 
made no attempt to exclude Type I AGN or QSOs from 
our analysis. In order to discover how significant an omis- 
sion this is, we have repeated our analysis excluding all 
th ose objects appearin g in the X-ray sele c ted c atalogs 
of lSzokoiv et~aLl (120041 ) and iTreister et al.l (|2008f ). This 
excludes 3.2 % (250/7840) of galaxies from our main sam- 
ple, and 6.5 % (96/1482) from our A/, > 10" M© sam- 
ple. Given these numbers, it is perhaps unsurprising that 



the exclusion of these objects does not greatly affect our 
results: 7tot and 7fi.ac change by —0.07 and -1-0.03, re- 
spectively. 

9.7. Quantifying Potential Systematic Errors — 
Summary 

How, and how much, do our results depend on our 
analytical methods and experimental design? We have 
now described a rather large number of sensitivity tests, 
designed both to identify which aspects of our experi- 
mental design are crucial in determining our results, as 
well as to estimate the size of lingering systematic errors 
in our results. The results of many of these tests are 
summarized in Table [H Clearly, these sorts of tests can 
provide a staggering array of 'metadata', offering insights 
to guide not only the interpretation of the present data 
and results, but also the design of future experiments. 

As we have shown in ^8A\ the systematic uncertainties 
in the 7s due to discrepancies between the treatment 
of z w and z ;2> galaxies is on the order of ±0.2; 
this is the single greatest potential source of systematic 
uncertainty. Since we have treated these two samples in 
a uniform manner, we do not consider this as a lingering 
source of systematic uncertainty. 

The dominant sources of systematic uncertainty are, in 
order of importance: the method of deriving photometric 
redshifts (A7tot ~ 0.15, A7frac ~ 0.15); photometric cali- 
bration errors (A7tot ~ 0.05; A7fi.ac ~ 0.10); incomplete- 
ness (A7tot ~ 0.1, A7fi-ac ~ 0); the choice of photometric 
redshift template set (A7tot ~ 0.06; A7fiac ~ 0.06); pho- 
tometric methods (A7tot ~ 0.07; A7frac ~ 0); and stellar 
mass estimates (A7tot ~ 0.05; A7fiac ~ 0.05). We can 
also, for example, dismiss incompleteness due to our K 
S:N criterion, contamination of the sample by QSOs, and 
minor details of our photometric redshift calculation as 
significant sources of systematic uncertainty. 

In summary, we estimate the systematic errors on the 
measured values to be A7tot = 0.21 and A7fi.ac = 0.20; 
these values have been obtained by adding the dominant 
sources of systematic error in quadrature. We note that 
previous studies have not generally taken (all) sources of 
systematic error into account in their analysis. 

10. A FINAL INDEPENDENT CONSISTENCY CHECK: 
QUANTIFYING THE EVOLUTION OF 
BRIGHT, RED GALAXIES WITHOUT REDSHIFTS 

Our goal in this paper has been to quantify the growth 
of massive galaxies in general, and of massive red galax- 
ies in particular, over the past 10 Gyr, or since z 2. We 
have now seen — not completely unexpectedly — that the 
use of photometric redshifts is potentially a major source 
of systematic error in such measurements. For this rea- 
son, in this section we present a complementary and com- 
pletely independent measurement which does not rely on 
redshift information at all. 

10.1. Selecting Red Galaxies in Redshift Intervals 

Figure [13] shows evolutionary tracks for a passively 
evolving stellar population, in terms of observed col- 
ors in the MUSYC E CDFS bands. Specifically, we 
have used Pegase V2.0 (|Le Borgne fc Rocca-Volmerangd 
|2002| ) models with a short star formation episode (e- 
folding timescale of 100 Myr), beginning with Z ~ 0.004 
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Fig. 13. — Evolutionary tracks, in terms of observed colors, for 
a passively evolving stellar p opulation — Each track is for a syn- 
thetic Pegase V2.0 IILe Borg ne & Rocca-Vol mcrang c 2002) stellar 
population that evolves passively from Zform = 5; each of these 
tracks describe an approximate 'red envelope' for the true 2^2 
color-redshift distribution. The sharp rise in each track is due to 
the Balmer and/or 4000 A breaks becoming redshifted between the 
two observed filters. The shaded regions show where each of the 
four observed color criteria given will select red galaxies ; taken to- 
gether, they therefore provide a means for selecting passive galaxies 
in rough redshift bins (see also i|10.2l l. 



TABLE 4 

Quantifying the evolution of bright, red galaxies 
without redshifts 



Color Criterion 


Redshift Prevalence 








2low 


(2) Observed 


Expected 


Relative 




(1) 


(2) 


(3) (4) 


(5) 


(6) 


B - 


R > 1.6 


0.1 


0.3 0.47 / □" 


0.76 / □" 


62 % 


V - 


/ > 2.0 


0.4 


0.6 0.80 / □" 


1.80 / □" 


45 % 


R- 


2 > 1.5 


0.6 


0.9 0.84 / □" 


1.94 / □" 


43 % 


I - 


J > 1.4 


1.1 


1.4 0.37 / □" 


1.29 / □" 


29 % 


J - 


K > 1.4 


2.0 


2.6 0.35 / □" 


1.55 / □" 


23 % 



Note. — We have used a tiered set of color selection criteria 
(Col. 1) to select red galaxies in approximate redshift intervals. 
Here, it is understood that redder criteria are given primacy; for 
example, objects counted in the second row may or may not satisfy 
the first criterion, but do not satisfy any of the three subsequent 
criteria. We give the minimum and mean redshifts of color-selected 
galaxies from the synthetic, 'passive evolution' catalog in Col.s (2) 
and (3), respectively. The observed prevalence of color-selected 
objects in the MUSYC ECDFS catalog are given in Col. (4); Col. 
(5) gives the same quantity for the synthetic catalog. The ratio of 
the two is given in Col. (6), which gives an estimate of the rela- 
tive number of (potentially) passively evolving galaxies at different 
redshifts. 



at Zform = 5, and ending with Z « Zq. This (approxi- 
mately) maximally old model describes a 'red envelope' 
for the colour-redshift relation for observed galaxies: at a 
given redshift only extreme dust extinction can produce 
observed colors redder than these tracks. 

As is evident in this figure, as the Balmer and/or 4000 
A breaks are first redshifted between the observed bands, 
there is a sharp rise in the observed colour — by se- 
lecting objects that are very red in a certain colour, 
it is thus possible to select galaxies that lie beyond a 
certai n redshift. This is completely analogous to the 
ERQ (lMcCarth v"2004'. and references therein) or DRG 
(jFranx et all 120031 selection criteria for red galaxies at 
moderate-to-high redshifts. The particular selections we 
have adopted in Figure [T3l (and given in Col. 1 of Table 
H]) translate to minimum redshifts of approximately 0.16, 
0.40, 0.66, 1.12, and 2.05. 

By applying several of these colour criteria in concert, 
with the reddest criteria given primacy, it is then pos- 
sible to select red galaxies in rough redshift intervals: 
the highest redshift galaxies are selected by the criterion 
{J~K) > 1.4; then, the next highest redshift galaxies are 
then selected from the remaining set {ie. (J — K) < 1.4) 
by the criterion (/ — J) > 1.4, and so on. In the current 
context, this can be thought of as a two-colour, binned 
photometric redshift. 

The prevalence of objects selected in this way are given 
in Col. 4 of Tabled To our hmit of < 22, we find 386, 
655, 690, 304, and 286 objects selected by this tiered set 
of colour criteria. 

10.2. Interpreting the Numbers of Colour Selected 
Galaxies 

It is clear that the exact redshift range over which a 
galaxy might satisfy a given color criterion depends on 
that galaxy's SED: as well as passively evolving galax- 



ies, these selections will also include galaxies whose red 
observed colors are due to, e.g., dust obscuration or con- 
siderably higher redshifts. Thus, the details of the evolv- 
ing, bivariate color-magnitude distribution is folded into 
the numbers of color-selected galaxies. For this reason, 
we are forced to interpret the numbers of color-selected 
galaxies with reference to a particular model for the evo- 
lution of red galaxies. 

This is done as follows: we construct a mock catalog 
in which galaxies' luminosities are distributed accord- 
ing to the z « luminosity function for red galaxies, 
which we have determined from the low-z sample, as an- 
alyzed in Appendix jAl This lumin osity function is nearly 
identical to that of iBlanton et al I ((2005cf) : our results do 
not change significantly a dopting th e red g ala xy luminos- 
ity fun ctions from either iBell et aP (|2003l ) or ICole et aTl 
(|2001h . Galaxies are placed randomly (ie. uniform co- 
moving density) in the volume z < 5. Next, we generate 
synthetic photometry fo r each object in the catalog using 
the s et of Pegase V2.0 ()Le Borgne fc Rocca-Volmerangd 
|2002|) models shown in Figure [T31 which are essentially 
passively evolving from Zform = 5. We also include typ- 
ical errors for the MUSYC ECDFS catalog, to approxi- 
mately account for photometric scatter. Finally, we ap- 
ply our colour selection criteria to this catalog, exactly as 
for the observed ECDFS catalog. The predicted preva- 
lence of colour selected galaxies, to be compared to those 
observed, are given in the Col. 5 of Tabled 

Whereas all red galaxies are assumed to evolve pas- 
sively in the synthetic reference catalog, the real color 
selected samples will include dusty or high-redshift star- 
forming galaxies: while all passive galaxies are red, not 
all red galaxies are passive. As we argue in §8.11 the 
number of color selected galaxies can therefore be used 
to place an upper limit on the number of passively evolv- 
ing galaxies, since these color criteria should select a 
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complete but contaminated sample of genuinely passive 
galaxies. 

The results presented in Table [3] thus suggest that the 
number density of bright, passively evolving galaxies has 
increased by a factor of at least ~ 2 since z ~ 1. More- 
over, we find that the number density of passive galax- 
ies in the range 1 < z < 2 is at most ^ 25 % of the 
present day value. Taken at face value, this would sug- 
gest that at most 1/4 present day red sequence galaxies 
can have evolved passively from z ~ 2; the remainder 
were still forming, whether through star formation or by 
the hierarchical assembly of undetected, and thus fainter, 
subunits, or some combination of the two. 

Finally, we emphasize the essential simplicity of this 
analysis, and thus its suitability for comparisons between 
different datasets, which may use different strategies for 
the computation of photometric redshifts, restframe col- 
ors, stellar masses, and so on. 

10.3. Tying it all Together 

Our final task is to compare the results of this section 
with those in iJH we do this in Figure [T4l In this plot, 
the abscissa is the effective wavelength of the bluer of the 
two filters used in each color criterion; the ordinate is the 
number of color selected galaxies observed in the MUSYC 
ECDFS catalog, relative to the expectation for passive 
evolution. The yellow squares refer to the Zform = 5 
model shown in Figure fT3l and presented in Table ID The 
blue triangles and red pentagons show how these results 
would vary had we assumed Zform = 4 and Zform = 6, 
respectively. 

The specific redshift range that is selected by each color 
criterion depends on the particular tracks used. For a 
given color criterion, the mean redshift of selected sources 
(derived from the synthetic 'passive evolution' catalog) 
thus varies for different choices of Zform: approximate 
transformations are given for each scenario at the top of 
the Figure. 

Finally, the solid line shows our fit results for the 
•Zphot ^5 2 number density evolution of red galaxies, 
shown in Figure [TO] . In order to put this curve on this 
plot, we have used the mean redshift of (synthetic, pas- 
sive) color-selected sources assuming passive evolution 
from Zform = 5. Given the very different assumptions 
and methods lying behind these two results, the agree- 
ment is very impressive. 

Both analyses clearly indicate that the number of red 
galaxies at z « 0.7 is approximately 50 % of the present 
day value; at z « 1.5, it is approximately 25 %. In other 
words, at most, approximately 1/2 local red sequence 
galaxies were already 'in place' by z « 0.7, and have 
evolved passively since that time; at z « 1.5, at most 
1/4 were in place. 

11. SUMMARY AND CONCLUSIONS 

We have presented the color-magnitude (Figure |4]) and 
color-stellar mass diagrams (Figure [5]) for z < 2, based 
on a sample of < 22 galaxies drawn from the MUSYC 
catalog of the ECDFS. On the basis of the ten band 
SEDs, we achieve a photometric redshift accuracy of 
fJz = 0.036 (Figure [T|); this figure represents the cur- 
rent state of the art for broadband photometric look- 
back surveys. We have used an empirical argument to 
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Fig. 14. — Quantifying the number evolution of bright, red 
galaxies without redshifts — The observed numbers of color-selected 
galaxies, relative to the number expected for passive evolution, 
plotted as a function of the effective wavelength of the bluer filter 
used in the selection criterion; the criteria themselves are given 
in Table |4] The ordinate in this plot is essentially the number of 
galaxies with colors that are consistent with passive evolution from 
^formi relative to the number of 2 red sequence galaxies; we 
have plotted results assuming Zform = 4 {blue triangles), 5 {yellow 
squares), and 6 {red pentagons). The mean redshift of color selected 
galaxies, as derived from the synthetic 'passive evolution' catalog, 
does depend mildly on the choice of evolutionary scenario. The 
mean redshift of galaxies selected by each criterion, (z), is given 
separately for each choice of 2form ^-t the top of the Figure. The 
minimum redshift is more robust; this is given as Zmin- Lastly, the 
solid line and shaded region show the fit to the Zpjujt ^ 2 evolution 
in the number density of massive red galaxies, derived in JS] (see 
Figure [To] l: the shaded region shows the Ict uncertainty in the fit. 
(To put this relation on this plot, we have used the Zform = 5 cat- 
alog to map each color criterion to {z).) The agreement between 
these two very different calculations is impressive. 



demonstrate that our main galaxy sample is approxi- 
mately complete (volume limited) for > 10^^ and 
2phot < 1-8 (Figure ED. 

Based on the joint color-stellar mass-redshift distribu- 
tion of this mass-selected sample, we make the following 
conclusions: 

1. The color distribution of the massive galaxy popu- 
lation is well described by the sum of two separate 
Gaussian distributions for Zphot ^1-2 (Figure [6]). 
Beyond this redshift, the depth of our NIR data 
makes it impossible to identify distinct red and blue 
subgroups from within the general massive galaxy 
population on the basis of the observed color distri- 
bution. The question as to the existence or other- 
wise of a Zphot ^1-5 red sequence will require data 
approximately an order of magnitude deeper than 
our own (Figure [7]). 

2. The colors of red sequence galaxies have become 
progressively redder by ~ 0.5 mag since Zphot ~ 1.1 
(Figured]). Making a linear fit to the observed evo- 
lution, we find A{u - r) cx (-0.44 ± 0.02) Zphot- 
Simple models can only describe the observations 
assuming sub-solar metallicities (Figure [9]) : it re- 
mains a challenge to consistently describe the ob- 
served colors and the mass-metallicity relation. 
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3. While the number density of massive galaxies 
evolves mildly for Zphot < 2, we see strong evo- 
lution in the red galaxy fraction. That is, the mas- 
sive galaxy population appears to be changing more 
than it is growing. We have quantified this evolu- 
tion using parametric fits of the form (1 + z)^ . For 
Zphot < 1.8, we find 7tot = -0.52 ± 0.f2(±0.20), 
7rcd = -1.60± 0.14(±0.2f), and 7ft.ac = -1.06 ± 
0.16(±0.2f). 

The systematic errors (given in brackets above) have 
been derived on the basis of a whole raft of sensitivity 
analyses, and are due primarily to photometric calibra- 
tion errors and systematic effects arising from the pho- 
tometric redshift calculation (f}9|). 

Finally, in iJTOl we showed that these results are com- 
pletely consistent with an independent analysis based 
only on directly observed quantities; that is, without de- 
riving redshifts, etc., for individual galaxies. 

The two major advances in this work are the quan- 
tification of the Zphot ^ 1.2 color evolution of the red 
galaxy sequence as a whole, and the quantification of the 
^phot ^ 2 evolution of the red fraction among > 10^^ 
Mq galaxies. 

Knowing that the vast majority of 10^^ Mq galaxies in 
the local universe fall on the red sequence (Figure ^ , we 
can identify our massive galaxy sample as the immediate 
progenitors of local red sequence galaxies. Further, by 
extrapolating the observed Zphot < 1.4 color evolution 



of the red sequence to higher redshifts, we can identify 
our 'red' galaxies as those which potentially have already 
found their place on the red sequence: the high redshift 
analogues of local red sequence galaxies. However, sim- 
ply selecting 'red' galaxies may include a significant num- 
ber of galaxies whose red colors are due to heavy dust 
obscuration, we argue that our red galaxy measurements 
provide an approximate upper limit on the number of 
passive or quiescent galaxies. 

With this assumption, our results suggest that at least 
1/2 of all Af, > 10"^^ Mq galaxies joined (or rejoined) 
the red sequence only after z ^ 1, and that at most 1/5 
massive galaxies have resided on the red sequence since 
z ~ 2. These results provide new constraints for quench- 
ing models , such as quasa r and "radio mo de" AGN feed- 
back (CrotoZetaDlSol; |CattanioIiE3 12006). Estab- 
lishing which of our 'red' galaxies are genuinely 'red and 
dead' offers the opportunity to considerably tighten these 
constraints. 
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and participants of the several workshops hosted by the 
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APPENDIX 

A. THE Z = COMPARISON POINT 

In this Appendix, we describe the process by which we have derived the z « number density of massive galaxies 
in general, and of red galaxies in particular, from the low-z subsample of the New York University (NYU) Value 
Added Galaxy Catalog (VAGC) fr om the Sloa n Digital Sky Survey (SDSS) data release 4 (DR4); this catalog has 
been compiled and described by .Blanton et al.l (|2005b) . We will refer to this sample simply as the 'low-z' sample, as 
opposed to the 'high-z', K-selected sample in the main text. The overarching concern is uniformity in the analysis 
of the low-z and high-z samples, where possible and appropriate, in order to make as fair as possible a comparison 
between the high-z and low-z data. Our discussion proceeds in three parts: first, we describe our analysis of low-z 
galaxies; we then describe our characterisation of the galaxy red sequence; finally, we give our derived number density 
of massive galaxies at z « 0, including mimicking the effects of the random photometric redshift errors typical for the 
z 3> sample. 

A.l. Low-z Galaxy Analysis 

For each galaxy in the low-z sample, we have constructed a ugriz SED using 'model' magnitudes (measured by fitting 
either an exponential or r^/"* profile in each band, with structural parameters derived from the i band image) as given 
in the low-z catalog. In analogue to high-z galaxies, these SEDs are then scaled up to the 'Petrosian' r magnitude. The 
Petrosian apertures are designed to measure a fixed fraction of a galaxy's light, irrespective of brightness or distance, 
but with some dependence on the light distribution: neglecting the effects of seeing, the Petrosian ape rture captures 
99 % of the total light for an exponential profile, compared to > 80 % for a de Vaucouleurs profile (jBlanton et al.l 
l2005bf) . We make no attempt to corre ct for missed flux, but note that similar effects are present at a similar level in 
our own data (see i j9.2.2p . Following (jBlanton et al.|[2b05bf ). we use the following factors to correct the basic SDSS 
calibration to AB magnitudes: (-0.042, -^0.036, +0.015, +0.013, -0.002) for {u, g, r, i, z). 

We have derived restframe photometry and stellar masses for low-z galaxies using exactly the same machinery 
as for the high-z sample, adopting the heliocentric redshifts g iven in the low-z catalog. Our interpolated restframe 
ugriz photometry agrees quite well with that given by Blanton et al.l ()2005bl ) in the low-z catalog. (Whereas we 
use galaxy color-color relations, based on six empirical galaxy template spectra, to interpolate restframe photometry, 
iBlanton et al.l (|2005b[ ) determine restframe fluxes by fitting a non-negative linear combination of five carefull y chosen 
synthetic spectra, and then integrating the best fit spectrum.) On average, in comparison to IBlanton et al.l ([^05b), 
our u — r colors are ^ 0.02 mag bluer for blue galaxies and ~ 0.03 mag redder for red galaxies. We calculate the 
distance modulus using the proper motion corrected redshifts given in the catalog. Finally, we have derived stellar 
masses based on the interpolated restframe B and V fluxes, using Equation [TJ 
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A. 2. The Red Galaxy Sequence at z — 

Figure ED shows our CMD (left) and CM*D {right) diagrams for the general z « field galaxy population; in each 
panel, the logarithmic greyscales show the data density. We note that these two plots are very nearly equivalent: the 
CM^D can be seen as simply a 'sheared' version of the CMD, with the bluer galaxies dragged further towards the left. 
A separate red galaxy sequence is easily discernible in each panel. 

We have characterized this red galaxy sequence, as a function of magnitude or of mass, by taking a narrow slices 
in either absolute magnitude or stellar mass, and fitting double Gaussians to the color distribution in each bin. The 
results of these fits are shown in each panel of Figure [AT] as the red and blue contours. This gives some indication of 
the relative numbers of red and blue galaxies as a function of magnitude/mass. In order to characterize the separate 
red and blue populations, we also show the center and width of each Gaussian fit as the yellow and cyan points and 
error bars. 

In order to determine the slope of the high hmrinosity/mass end of the red sequence OMR and CM*R, we have 
simply made linear fits to the centers of the Air < —21 and Af* > 10^°'^ Mq 'red' distributions. The results of these 
fits are: 

{u-r) = 2.612 - 0.090 x {Mr + 22) = 2.573 + 0.237 x logio(M*/10" Mq) . (Al) 

These relations are shown as the solid black lines in each panel. Finally, in order to make some distinction between 
'red' and 'blue' galaxies which could be easily applied to the high-z sample, we simply use a cut 0.25 mag bluer than 
this relation. In each panel, this cut can be seen to approximately coincide with the so-called 'green valley' between 
the red and blue galaxy populations. 

Finally, in the inset of each diagram, we show the color distribution for the brightest {Air < ~21) or most ma ssive 
(M* > 10^^) galaxies, after subtracting away the color-magnitude or color-stellar mass relations given by Equation lAll 
(See also El) We also show double Gaussian fits to these distributions, which are the basis for the z ^ points shown 
in Figures [H and [H We note in passing that, at least for —19 > Mr > —21.5, the width of the color distributions of 
red and blue galaxies is not a strong function of magnitude: A{u — r) ^ 0.12, and A(m — r) w 0.20 for red and blue 
galaxies, respectively. 

A. 3. The Number Density of Massive Galaxies at z « 

For our purposes (see ^JH]), the crucial final quantity is the number density of massive galaxies at z « 0. We find that 
the number density of > 10^^ Mq galaxies is (6.64 ±0.38) x 10"'* Mpc""^; for (u, r) red galaxies, the number is 
(6.07 ± 0.36) X 10"'' Mpc~'^. This measurement is not significa ntly affected by e ither volume or s urface brightness 
incompleteness. Compared to the more sophisticated analyses of iBell et al.l (|2003[ ) and ICole et al.l ()2001[ ). our value 
for the number density of massive galaxies is 10 % higher and 25 % lower, respectively. 

As we have repeatedly stressed, our prime concern is uniformity in the analysis of the low- and high-z samples. 
Photometric redshift errors are a major source of uncertainty for the high-z sample — how would comparable errors 
affect the lowz measurements? Given the redshift range of our low-z galaxies (0.003 < z < 0.05) and the estimated 
photometric redshift errors among high-z galaxies (Az/(1 + z) ^ 0.035), it would be inappropriate to simply apply 
typical redshift errors and repeat our calculations: this would effectively throw away all distance information. Instead, 




-18 -19 -20 -21 -22 -23 9.5 10.0 10.5 11.0 11.5 12.0 

absolute magnitude, M,. stellar mass, logi(,(M,/Msoi) 

Fig. Al. — The 'low z' sample from the NYU VAGC of the SDSS (DR4), used to derive the z ^ comparison point for the main, 'high z' 
sample discussed in the text — Color— magnitude {left) and color— stellar mass (right) diagrams — The greyscale shows the density of points. 
The errors bars show the center and width of the color distributions of the red and blue populations, obtained from double Gaussian fits in 
narrow magnitude/mass slices. These fits to the populations are also represented by the logarithmic (0.5 dex) contours, to give some idea 
of the data density of red/blue galaxies. 
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what we have done is to apply the typical uncertainties on Mr, {u — r) and Af* due to photometric redshift errors, as 
shown in Figured For simplicity, we have ignored correlations between these errors. 

In line with our earlier findings, including the effects of photometric redshift errors makes less of a difference to 
measurements based on stellar masses than to those based on absolute magnitudes. Our fits to the red sequence 
become: 

{u-r) = 2.608 - 0.113 x {Mr + 22) ^ 2.577 + 0.223 x logio(Af,/10" Mq) . (A2) 

Our measurements of the z « number density of galaxies more massive than 10^^ M© becomes (6.93 ±0.38) x 10^** 
Mpc-3 in total, and (5.86 ±0.35) x 10"'' h?' Mpc'^ for red sequence galaxies. That is, the effects of the redshift errors 
expected among z ;» galaxies affect the measurement of the z sa number density by less than 5 %. 

These numbers, approximately accounting for the Eddington bias due to photometric redshift errors (as they apply 
to the z ^ sample), are then what have used as a local reference value to compare to the z data. 

B. A DETAILED COMPARISON WITH COMBO-17 

In the ECDFS, the COMBO-17 and MUSYC datasets represent parallel analytical paths describing the same physical 
reality — indeed the COMBO-17 broadband observations form the foundation of the MUSYC raw optical data. Ideally, 
then, the two surveys' results should agree identically; instead, in comparison to COMBO-17, the MUSYC results 
suggest 40 % more massive (Af, > 10^^ M©) galaxies at 0.2 < Zphot < 0.8 — a degree of difference more on a par with 
that expected from field-to-field variation. In this Appendix, we identify the main cause of this discrepancy. 

B.l. Photometric Redshift Accuracy: MUSYC photometry with COMBO-17 redshifts 

The essential differences between the COMBO-17 and MUSYC datasets are COMBO-17's twelve medium bands, 
allowing much more precise photometric redshift determinations for z < 1, and the MUSYC z'JHK data, which open 
the door to the z > 1 universe. The first and most obvious concern is thus that the difference between the two surveys' 
•Zphot ^ 1 results are a product of their different photometric redshift accuracies. To test this, we have tried simply 
adopting the COMBO-17 photometric redshifts in place of our own determinations, and repeated our analysis. 

The results of this test are shown in Figure [BT] the main panel shows our primary result {viz., the number density of 
Af* > 10"'^^ Mq galaxies as a function of redshift; see Figure fTO|). The solid histograms show this trial re-analysis of the 
MUSYC data adopting COMBO-17 redshifts; this should be compared to the fiducial results from MUSYC {dashed 
histogram) a nd COMBO-17 {d otted histograms). We have derived these 'COMBO-17' results fro m the p ublic catalog 
presented bv lWolf et al.l (|2003[ ). supplemented with the stellar mass determinations used bv Borch et al.l (|2006l . 

At least for Zphot < 0.8, the use of COMBO-17 redshifts in place of our own does not have a great effect on our 
results: the solid histograms lie quite close to the dashed ones. Quantitatively, adopting COMBO-17 redshifts in place 
of our own leaves the number of Af* > 10"'^^ Mq galaxies over 0.2 < Zphot < 0.8 unchanged to within 1 %; the effect 
is slightly greater for the red population, leading to a red fraction that is 17 % higher over the same interval. This 
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Fig. B1. — Reanalysis of the MUSYC data adopting redshifts from COMBO-17 — Main panel: The measured number density of M, > 10^^ 
Mq galaxies for z ^1 based on the MUSYC photometry, but adopting redshifts from the COMBO-17 catalog in place of our own {solid 
histograms). These results should be compared to the fiducial results from MUSYC {dashed histograms) and from COMBO-17 {dotted 
histograms). The red and black histograms are for red and all galaxies, respectively. Where the COMBO-i7 and MUSYC results differ by 
40 %, the effect of using of COMBO-17 redshifts is less than 10 %; the difference between the two surveys' results cannot be explained by 
differences in the photometric redshifts. Left panels: the difference between the {from top to bottom) photometric redshifts, absolute B and 
V magnitudes, restframe {B — V) colors, and stellar masses used to produce the results shown in the main panel, in comparison to those 
found in the COMBO-17 catalog, plotted as a function of the COMBO-17 redshift. In comparison to COMBO-17, even using the same 
redshifts, the MUSYC data and/or data analysis make galaxies appear brighter and redder, and so more massive. 
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implies that, at least for Zphot ^ 1, our results are not seriously affected by our lower photometric redshift accuracy 
{a, « 0.03 for MUSYC, versus 0.02 for COMBO-17, for i? < 24 and Zspcc < 1). 

To the left of the main panel in Figure IbTI we show the difference in {from top to bottom) the photometric redshifts, 
absolute B, and V magnitudes, restframe {B — V) colors, and stellar masses used to produce the solid histograms, and 
those from the COMBO-17 catalog. In all cases, the '5' is in the sense of 'MUSYC re-analysis' minus 'COMBO-17 
catalog', and is plotted as a function of the COMBO-17 catalog redshift. Within each panel, we also give the biweight 
mean offset for each quantity, evaluated for sources with i? < 23 (to reduce random scatter due to photometric 
uncertainties) and 0.2 < Zphot < 0.5 (where the COMBO-17 photometry still samples restframe V) in the COMBO-17 
catalog. 

Looking now briefly at these panels, we see that for Zphot ^ 0.5, where the COMBO-17 value for the restframe 
V magnitude comes from an extrapolation of the best-fit template spectrum, there is a progressive offset between 
the restframe V luminosities inferred by the MUSYC and COMBO-17 data and analyses — even while using the same 
redshifts. Even for Zphot < 0.5, however, we see a systematic offset of 0.05 mag in {B — V) color. Finally, we note that 
the greater number of 0.6 < Xphot < 0.8 galaxies in comparison to 0.2 < Zphot < 0.6 that we see in the MUSYC data 
is also present in the COMBO-17 photometric redshift distribution; the difference is that these galaxies do not have 
M, > 10" Mq in the COMBO-17 catalog. 

We therefore conclude that the difference between the COMBO-17 and MUSYC results is not a product of the two 
surveys' different photometric redshift accuracies: the use of COMBO-17 rather than MUSYC redshifts affects our 
results only by a few percent. This leaves the possibilities that the difference between the two results is due to the 
different methods used to infer galaxies' restframe properties, or to differences in the data themselves. 

B.2. Derivation of Restframe Properties: MUSYC analysis of the COMBO-17 photometry 

Are the discrepancies between the COMBO-17 and MUSYC results due to different systematic effects inherent in 
our different methods for deriving restframe photometry and stellar masses? In order to address this concern, we have 
reanalyzed the COMBO-17 data using the procedures described in the main text. For this test, in order to isolate 
the effect of the different analytical methods, we also adopt the COMBO-17 redshift determinations; the difference 
between this test and the previous one is thus the use of the COMBO-17, rather than the MUSYC, photometry. 

The results of this re-analysis are shown as the solid histograms in Figure lB2l all other symbols and their meanings 
are as in Figure [BT] The results of the MUSYC re-analysis of the COMBO-17 data agree very well with the COMBO- 
17 analysis proper: the solid black histogram lies very near the dotted black histogram. Quantitatively, the MUSYC 
analysis of the COMBO-17 data leads to a red fraction which is 33 % higher than for COMBO-17's own analysis; the 
total number density of Af* > 10^^ Mq galaxies agree to better than 1 %. 

Comparing the MUSYC- and COMBO- 17-derived restframe properties — again, based on the same photometry and 
redshifts — we do see systematic differences. We make three specific observations: 1. there appear to be discrete redshift 
regimes where the MUSYC- and COMBO-17-derived restframe fluxes compare differently; 2. for 0.2 < Zphot S 0.35, we 
see a 'bimodal' offset in the MUSYC- and COMBO- 17-derived {B — V) colors, corresponding to red and blue galaxies; 
3. the MUSYC-derived restframe fluxes appear to be systematically brighter than those derived by COMBO-17, even 
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Fig. B2. — MUSYC analysis of the COMBO-17 photometry — The soHd histograms in the main panel show the measured number density 
of M* > 10^^ galaxies, applying the MUSYC analysis to the COMBO-17 photometry. In order to isolate the effect on our measurements due 
to the different methods for deriving rest frame properties, we continue to adopt the COMBO-17 redshifts for this test. All other symbols 
and their meanings are as in Figure IBll Applying the MUSYC analysis to the COMBO-17 photometry, we reproduce the COMBO-17 
result to better than 1 %; the difference between the COMBO-17 and MUSYC results cannot be explained solely by our different analytical 
methods. 
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for the same photometry. Remarkably, even despite these differences, our stellar mass estimates agree to within 0.04 
dex, with a scatter of just 0.11 dex. Finally, we note that the progressive offset in extrapolated V luminosities for 
•s^phot ^ 0.5 seen in the previous test is not seen here; that is, applied to the same data, the two techniques for 
extrapolating restframe photometry yield similar results. 

We therefore conclude that the difference between the COMBO-17 and MUSYC results in the ECDFS cannot be 
explained by differences in the analytical methods employed by each team: the MUSYC reanalysis of the COMBO-17 
data agrees with the COMBO-17 fiducial results. Instead, it seems that the different results are really in the data 
themselves. 

B.3. Photometric Calibration: MUSYC analysis of the recalibrated COMBO-17 photometry 

A direct, object-by-object comparison of COMBO-17 and MUSYC photometry reveals significant differences in the 
photometric calibration of the two surveys (Paper I). Specifically, we see differences of 0.00, 0.06, 0.08, 0.08, and 0.14 
mag between the MUSYC and COMBO-17 U BV RI-ha,nd photometry for galaxies, such that the MUSYC photometry 
is systematically brighter and redder. This has subsequ ently been confirmed to be due to an error in the photometric 
calibration of the COMBO-17 data (|Wolf et al.l l2008l ) . Can this difference in photometric calibration explain the 
diflterent results found by COMBO-17 and MUSYC? 

To address this issue, we have simply scaled the COMBO-17 photometry to match the MUSYC photometry, and 
repeated our analysis. The COMBO-17 team has calibrated each of the medium bands relative to the nearest broad 
band; accordingly we have scaled each of the medium bands by the MUSYC-minus-COMBO-17 offset for the nearest 
broad band. The results of this test are shown in Figure [B3l The MUSYC re-analysis of the recalibrated COMBO-17 
photometry agrees extremely well with the fiducial MUSYC results: the measured number density of M* > 10^^ Mq 
galaxies at 0.2 < Zphot < 0.8 agree to better than 1 %. 

There are two separate aspects to this re-calibration: galaxies are both brighter and redder in the MUSYC catalog 
than in COMBO-17. Since COMBO-17 measures total fiuxes from their R band image, the 0.08 mag offset between 
the COMBO-17 and MUSYC zeropoints makes galaxies appear 7.6 % brighter (and thus more massive) in the MUSYC 
catalog; this effect is responsible for approximately 70 % of the change in the measured number density. At the same 
time, the reddening of the SED shape implies a higher mass-to-light ratio after recalibration; this effect is responsible 
for the other 30 % of the change in the measurements. 

We therefore conclude that the primary cause for the disagreement between the results of the MUSYC and COMBO- 
17 surveys in the ECDFS is the different photometric calibrations of the two surveys: reanalyzing the COMBO-17 
data set, recalibrated to match the MUSYC photometry, the results agree with the MUSYC fiducial analysis. 
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Fig. B3. — Recalibration of the COMBO-17 photometry — the solid histograms in the main panel show the measured number density 
of Mt > lO^'^ galaxies, applying the MUSYC analysis to the COMBO-17 photometry, recalibrated to match the MUSYC photometry. 
Specifically, the UBVRI bands have been scaled up by 0.00, 0.06, 0.08, 0.08, and 0.14 mag; each of the medium bands has been scaled to 
match the nearest broad band. We continue to adopt the COMBO-17 redshifts for this test. All other symbols and their meanings are as in 
Figure [BT] The effect of this rescaling is to increase the measured values by a further 30 %, in comparison to the previous test. Together, 
the combined effect of this recalibration is nearly 50 %, almost completely explaining the difference between the COMBO-17 and MUSYC 
results. 



